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Abstract 

The convergence of the Bohzmann equaiton to the compressible Euler equations when the Knudsen 
number tends to zero has been a long standing open problem in the kinetic theory. In the setting 
of Riemann solution that contains the generic superposition of shock, rarefaction wave and contact 
discontinuity to the Euler equations, we succeed in justifying this limit by introducing hyperbolic 
waves with different solution backgrounds to capture the extra masses carried by the hyperbolic 
approximation of the rarefaction wave and the diffusion approximation of contact discontinuity. 
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As the fundamental equation in statistical mechanics, the Boltzmann equation takes the form of 
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where f{t,X,^) is the density distribution function of particles at time t with location X and velocity 5. 
In this equation, the physical parameter e > called Knudsen number is proportional to the mean free 
path of the interacting particles. 

It was known since its derivation that the Boltzmann equation is closely related to the systems of fluid 
dynamics, in particular, the systems of Euler and Navier-Stokes equations. In fact, the first derivation of 
the fluid dynamical components and systems from the kinetic equations can be traced back to the dates 
of Maxwell and Boltzmann. Their early derivations rest on some arguments as how the various terms in 
a kinetic equation balance each other. These balance arguments seem arbitrary to some extent. For this, 
Hilbert proposed a systematic expansion in 1912, and Enskog and Chapman independently proposed 
another expansion in 1916 and 1917 respectively. 

Either the Hilbert expansion or Chapman-Enskog expansion yields the compressible Euler equations 
in the leading order with respect to the Knudsen number e, and the compressible Navier-Stokes equations, 
Burnett equations in the subsequent orders. To justify these formal approximations in rigorous math- 
ematics, that is, hydrodynamic limits, has been proved to be extremely challenging and most remains 
open, in part because the basic well-posedness and regularity questions are still mostly unsolved for these 
fluid equations. 

The justification of the fluid limits of the Boltzmann equation is also related to the Hilbert's sixth 
problem, " Mathematical treatment of the axioms of physics" , in which it says that " The investigations 
on the foundations of geometry suggest the problem: To treat in the same manner, by means of axioms, 
those physical sciences in which mathematics plays an important part; in the first rank are the theory 

of probabilities and mechanics Thus Boltzmann's work on the principles of mechanics suggests the 

problem of developing mathematically the limiting processes, there merely indicated, which lead from 

the atomistic view to the laws of motion of continua Further, the mathematician has the duty to 

test exactly in each instance whether the new axioms are compatible with the previous ones." 

The goal of this paper is to justify the limiting process of the Boltzmann equation to the system 
of the compressible Euler equations in the setting of Riemann solutions. The Riemann problem was 
first formulated and studied by Riemann in 1860s when he studied the one space dimensional isentropic 
gas dynamics with initial data being two constant states. The solution to this problem turns out to be 
fundamental in the theory of hyperbolic conservation laws because it not only captures the local and 
global behavior of solutions, but also fully represents the effect of the nonlinearity in the structure of the 
solutions. It is now well known that for the system of Euler equations, there are three basic wave patterns, 
that is, shock wave, rarefaction wave and contact discontinuity. These three types of waves have essential 
differences, that is, shock is compressive, rarefaction is expansive, and contact discontinuity has some 
diffusive structure. Therefore, how to study the hydrodynamic limit of Boltzmann equation for the full 
Riemann solution that consists of the superposition of these three typical waves is still very challenging 
in mathematics. 

In this paper, by introducing two types of hyperbolic waves that carry the extra masses in the 
hyperbolic approximation of the rarefaction wave profile and the diffusive approximation of the contact 
discontinuity, we succeed in proving rigorously that there exists a family of solutions to the Boltzmann 
equation that converges to a Maxwellian determined by a Riemann solution consisting of three basic wave 
patterns when the Knudsen number tends to zero. Furthermore, a convergence rate is obtained in term 
of the Knudsen number. 

By coping with the essential properties of individual wave pattern, the hydrodynamic limit for a 
single wave was justified in the previous works separately. More precisely, by using the compressibility of 
shock wave profile, Yu [41] showed that when the solution of the Euler equations ()1.2p contains only non- 
interacting shocks, there exists a sequence of solutions to the Boltzmann equation that converge to a local 
Maxwellian defined by the solution of the Euler equations (II. 2[) uniformly away from the shock in any fixed 
time interval. In this work, a generalized Hilbert expansion was introduced, and the analytic technique 
of matching the inner and outer expansions developed by Goodman- Xin [16] for conservation laws was 
used. On the other hand, by using the time decay properties of the rarefaction wave, similar problem 
was studied by Xin-Zeng [3D]. Moreover, by using the diffusive structure in the contact discontinuity 
as for the Navier-Stokes equations, the hydrodynamic limit to the contact discontinuity was proved by 
Huang- Wang- Yang in [T9] . 

However, up to now, how to deal with the general Riemann solution that consists of all three basic 
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waves is still a challenging open problem. This is mainly due to the difficulty in handling the wave 
interactions and also unifying the different approaches in the analysis used for each single wave pattern. 
In order to overcome these difficulties in justifying the limit, our main idea in this paper is to introduce 
two families of hyperbolic waves, called hyperbolic wave I and II, that capture the propagation of the 
extra mass created by the approximate hyperbolic rarefaction wave profile in the viscous setting and the 
diffusion approximation of contact discontinuity. 

We now briefly explain why the two families of hyperbolic waves we introduced are essential for 
the proof. As in the previous works on the rarefaction wave in the setting of either Navier-Stokes 
equations or the Boltzmann equation, the approximate rarefaction wave is constructed as a hyperbolic 
wave profile. Therefore, we need to precisely capture the error in the second order of the approximation 
for the Boltzmann equation in term of the Knudsen number, that is, in the Navier-Stokes level. And this 
reduces to study the propagation of the extra mass induced by the viscosity and heat conductivity. For 
this, we introduce the hyperbolic wave I as a solution to the linearized system around the approximate 
rarefaction wave profile with source terms given by the viscosity and heat conductivity induced by the 
rarefaction wave profile to recover the viscous terms. We can show that the hyperbolic wave I decays like 
the first-order derivative of the rarefaction wave profile so that the decay properties given in Lemma [2?2] 
are good enough to carry out the analysis. 

The main difficulty comes from the approximation of the contact discontinuity. First of all, such an 
approximation, that is, 2-viscous contact wave, behaves like a diffusion wave profile as for the Navier- 
Stokes equations. Due to the lack of sufficient decay in e and the non-conservative error terms when taking 
the anti-derivative of the perturbation, we need to remove the leading error terms and non-conservative 
terms in such approximation before taking the anti-derivative. The hyperbolic wave II is constructed to 
remove these error terms due to the viscous contact wave approximation. Note that the construction 
of the hyperbolic wave II can not be done simply around the 2-viscous contact wave as the hyperbolic 
wave I for the rarefaction wave profile. Otherwise, the wave interaction terms thus induced will lead to 
insufficiently decay in e due to 2-viscous contact wave and it seems that these terms are essential in wave 
interactions. Instead, it is constructed around the superposition of the approximate 1-rarefaction wave, 
the hyperbolic wave I, the 2-viscous contact wave and the 3-shock profile as a whole. Thus some wave 
interaction terms can be absorbed in the hyperbolic wave II and the other wave interaction terms can be 
handled by some subtle and careful calculations. Due to the non-conservative terms and insufficient decay 
rates of e of error terms induced by the 2-contact wave, we can not use the anti-derivative technique to 
analyze the hyperbolic wave II. Since the derivative of 3-shock profile is negative and tends to infinity as 
the Knudsen number £ — > 0+, we have to impose the condition at the time t = T (see (|2.74|) below) for 
the linearized hyperbolic system of hyperbolic wave II so that the monotonicity of 3-shock wave is fully 
utilized. This idea is different from the previous stability analysis on the shock profile which is based on 
the anti-derivative technique. 

With the help of these two hyperbolic waves and the corresponding new estimates, we can justify the 
limiting process from the Boltzmann equation to compressible Euler equations for the generic Riemann 
problems by elaborate analysis after a hyperbolic scaling. 

We now formulate the problem. Consider the Boltzmann equation with slab symmetry 

ft + ^if.^^QUJ), (1.1) 

where ^ = (^17^2,^3) ^ ^-^^d x E R^. Here, the collision operator takes the form of 

Q{f,gm {fiOgiO + fiQgia - /(05(?*) - fi^*)9iO)Bi\i-Ud) d^*dn, 

where £,',C* are the velocities after an elastic collision of two particles with velocities before the 
collision. Here, 9 is the angle between the relative velocity £, — £,* and the unit vector 51 in = {SI G 
S'^ : (C ^ C*) • > 0}. The conservations of momentum and energy yield the following relations between 
the velocities before and after collision: 
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We will concentrate on the hard sphere model where the cross-section takes the form of 

s(ie-e*i,^) = i(e-e*,f^)i = ie-e*|cos^. 

On the other hand, it is noted that the analysis can be applied to at least hard potential if we can assume 
Lemma [2731 on the shock wave profile holds true. 

As we mentioned earlier, formally, when the Knudsen number e tends to zero, the limit of the Boltz- 
mann equation (jl.ll) is the system of compressible Euler equations that consists of conservations of mass, 
momentum and energy: 



pt + {pui)x = 0, 

{pui)t + {pu\ +p)x = 0, 

{pui)t + {puiUi)x =0, i = 2,3, 

[p{e+^)]t + [pu,iE + l^)+pu,] 



(1.2) 



where 



p{t,x) 
pui{t, x) 



p(e + ^)(i,:r)= / 



f^iOfit,x,Od^, J -1,2,3, 



(1.3) 



Here, p is the density, u = (mi,M2,M3) is the macroscopic velocity, e is the internal energy, and p — RpO 
with R being the gas constant is the pressure. Note that the temperature is related to the internal 
energy by e = and (pi{£_){i = 0, 1, 2, 3, 4) are the collision invariants given by 



'^0(0 = 1, ¥^^(0 (^ = 1,2,3), ^4(0 =2 1^1'' 



(1.4) 



that satisfy 



R3 



ip^{OQi9l,92)d^^0, for z = 0,1, 2, 3, 4. 



Instead of using either Hilbert expansion or Chapman-Enskog expansion, we will apply the macro- 
micro decomposition introduced in [28]. For a solution /(i,x,^) of (jl.ip . set 

f{t,x,O^M{t,x,0 + G{t,x,0, 

where the local Maxwellian M(i, x,£,) — M[p „ (^) represents the macroscopic component of the solution 
defined by the five conserved quantities, i.e., the mass density p(t,x), the momentum pu(t,x), and the 



total energy p{e 



^'){t,x) given in (|1.3p . through 

M = M[p,„,e](t,a;,0 



p{t,x) 



(1.5) 



y^{2nRe{t,x)f 

And G(t,x,^) represents the microscopic component. 

From now on, the inner product of gi and 172 in i|(R-'^) with respect to a given Maxwellian M is 
denoted by: 



(51,32) 



M 



^giiMOd^. 

R3 M 



(1.6) 



If M is the local Maxwellian M defined in (jl.Sp . the macroscopic space is spanned by the following five 
pairwise orthogonal base. 



fp 

X^(0 = ^^^-^^^ for i = 1,2,3, 
X4(0^^(^^^3)M, 



(1.7) 



RB 

- (x»,Xj>=%, i, J = 0, 1,2,3,4. 
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For brevity, if M is the local Maxwellian M, we will simply use (•, •) to denote (•, •)m- By using the above 
base, the macroscopic projection Pq and microscopic projection Pi can be defined as 

4 

Note that a function g{^) is called microscopic if 

mv^iOd^^0, z = 0,1, 2, 3, 4, 



where again ipi(^) represents the collision invariants. 

Notice that the solution f{t,x,£,) to the Boltzmann equation (jl.ip satisfies 

Po/ = M, Pi/ = G, 

and the Boltzmann equation (jl.ll) becomes 

(M + G)t+ei(M + G), = i[2Q(M,G) + g(G,G)]. 



(1.8) 



By integrating the product of the equation (|1.8I) and the collision invariants (pi{^){i — 0, 1,2,3,4) with 
respect to ^ over R'^, one has the following system for the fluid variables (p, u, 9): 



pt + ipui)x = 0, 
{pui)t + {puf +p).. 



{pui)t + {pUlUi)^ = - 



iiUG^dS,, 1 = 2,3, 

2 



(1.9) 



Note that the above fluid-type system is not self-contained and one more equation for the microscopic 
component G is needed and it can be obtained by applying the projection operator Pi to (jl.Sp : 

G, + Pi(eiM,) + Pi(eiG,) = i [LmG + Q(G, G)] . (1.10) 

Here Lm is the linearized collision operator of Q{f, /) with respect to the local Maxwellian M given by 

Lm5 = 2g(M,5) = g(M,.9) + Q{g,m)- 
Note that the null space OT of Lm is spanned by the macroscopic variables: 

X, (e), J =0,1, 2, 3, 4. 

Furthermore, there exists a positive constant 5 > such that for any function g{£^) e ^Jl^, cf. [T7] . 

(5,Lm5) < -'^{v{\^\)g-,g)-, 

where v{\£\) = 0(1)(1 -I- |^|) is the collision frequency for the hard sphere model. 

Consequently, the linearized collision operator Lm is a dissipative operator on L^(R'^), and its inverse 
Lj^^ is a bounded operator on DT^. It follows from (|1.10p that 



with 



G = eLMi[Pi(eiM,)]+n, 



n = Lj^i[e(Gt + Pi(aG,)) - Q(G, G) 



(1.11) 
(1.12) 
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Plugging (|l.lip into (|1.9|) gives 



pt + {pui)^ = 0, 

{pui)t + {puiUi)^ = £{p{d)uix)x - J (.iCiUxdS^, 1 = 2,3, 

3 



(1.13) 



where the viscosity coefficient fj,{9) > and the heat conductivity coefhcient k{6) > are smooth functions 
of the temperature 6. Here, we normahze the gas constant i? to be | so that e = 9 and p — ^pO. 

Since the problem considered in this paper is one dimensional in the space variable a; G R, in the 
macroscopic level, it is more convenient to rewrite the equation (jl.ip and the system (|1.2p in the La- 
grangian coordinates. For this, set the coordinate transformation: 



{t,x)^ [t, / p{T,y)dy ~ {pui){T,y)dT], 

^ '(0,0) ^ 



(1.14) 



where jj^ fdy + gdr represents a line integraton from point A to point B on R+ x R. Here, the value of 
the integration is unique because of the conservation of mass. 

We will still denote the Lagrangian coordinates by {t, x) for the simplicity of notations. Then (|l.ip 
and (|1.2p in the Lagrangian coordinates become, respectively. 



ft 



fx H~ fx 
V V 



-Q{f, /), 



(1.15) 



and 



( vt- uix = 0, 

Ult +Px =0, 

Uit = 0, i = 2, 3, 



(1.16) 



Moreover, ([0|) - (ITTn)) take the form of 

Vt - Ulx = 0, 

Ult + Px ^ - J £,fGxd^, 
Uzt = - I ^i^iGxd^, i = 2,3, 



^ {0+\±)^ + ipui)x = -Jl^i\eG. 



d^, 



(1.17) 



with 



Gt - — + -Pi(^iM,) + -Fii^iGx) = ^(LmG + Q(G, G)), 

V V V e 



G = eLj^i(ipi(aM,)) + Hi, 



Hi = L^[e{Gt - ^Gx + iPi(eiG,)) - Q(G, G) 



(1.18) 

(1.19) 
(1.20) 
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and 



Vt - Ulj; = 0, 

uit +Px = — ( uixjx - / S,iUixd^, 

3 V J 



ili^^lxdi, i = 2,3, 



(1.21) 



2 i; 3 w 

i=2 ^ 

The Riemann problem for the Euler system (|1.16p is an initial value problem with initial data 

I a\u n ^ / (^-''"-'^-)' < 0, 
^ ' ' ' ^ \ (w+,u+,6'+), x>0, 

where, m = (ui, 1*2,^3), u± = (mi±,0, 0) and w± > 0, ui±,0± > are constants. It is known that the 
generic solution to the Riemann problem consists of three waves that propagates at different speeds, 
that is, shock, rarefaction wave and contact discontinuity, cf. [lOj [27]. We denote this solution by 
(F,f7,e)(t,x). Note that tj = (f7i,0,0). 

Given the right end state the following wave curves for the left end state {v,ui,9) in 

the phase space are defined with w > and 9 > {) for the Euler equations (jl.l6p . 

• Contact discontinuity curve: 



CD{v+,ui+,e+) = {{v,ui,e)\ui = = p+,w ^ v+). 

«-Rarefaction wave curve {i — 1,3): 



i?i(u+,wi+,6'+) := I {v,ui,( 



V < V+, Ul = Wl + 



X^{r],s+)dr], s{v,9) = 



(1.22) 



(1.23) 



where s+ — s(w+,6'+) and Xi — Xi{v,s) is the i-th characteristic speed of (|1.16p . 
• i-Shock wave curve (i — 1,3): 



-Si(w+ - u) - (ui+ - Wl) = 0, 

'Si{ui+ - Ul) + (p+ - p) = 0, and Xi+ < Si < Aj_ 
-Si{E+ - E) + {p+ui+ - pui) ^ 0, 



(1.24) 



l"+r 



204 
3v^ 



Xi± = Xi{v±,9±) and Si is the i— shock speed. 
For definiteness, we consider the case when the solution to the Riemann problem is a superposition 
of a 1-rarefaction and a 3-shock wave with a contact discontinuity in between, that is, (t;_,ui_,^?_) G 
Ri-CD-S3{v+,ui+,9+). Then there exist uniquely two intermediate states (w*, ui*, 0*) and (w*,u|,6'*) 
such that {v-,ui-,9-) S i?i(w*,ui*,0*), (v*,ui*,9*) G CD{v* ,u\,9*) and (i;*,u|,^*) e S'3(u+, 0+). 
Hence, the wave pattern {V , U,£){t,x) can be written as 




1 1 1 



(t,a;) 



1* 



, f/, = 0,(1 = 2,3), 



(1.25) 



E^ +E* 



where {v^'^ ,u\^ ,9'^'^){t,x) is the 1-rarefaction wave defined in (|1.23l) with the right state {v+,ui+,9+) 
given by {v^,,Ul^,,9^,), {v'^'^,u1'^,9'^'^){t,x) is the contact discontinuity defined in (|1.22p with the states 
{v-,ui-,9-) and {v+,ui+,9+) given by (w*,wi*,6'*) and {v*,ul,9*) respectively, and {v^^ ,u'^i'^ ,9''^){t,x) 
is the 3-shock wave defined in ()1.24p with the left state (u_,wi_,6'_) given by {v* ,ul,9*). 
Consequently, we can define 



e(t,x) = (£(i,x)-^%^). 



(1.26) 



7 



Due to the singularity of the rarefaction wave at i = 0, in this paper, we consider the problem in the 
time interval [h,T] for any small fixed /i > up to any arbitrarily fixed time T > 0. To investigate the 
interaction between the waves and the initial layer is another interesting topic that will not be discussed 
here. With the above preparation, the main result can be stated as follows. 

Theorem 1.1 Let {V ,Q){t, x) be a Riemann solution to the Euler equations which is a superposition 
of a 1 -rarefaction wave, a 2- contact discontinuity and a 3- shock wave, and 6 = \{v^ — V-,U-^-—u^,9+ — 0-)\ 
be the wave strength. There exist a small positive constant Sq, and a global Maxwellian = M[^^ u^ ^^j 
such that if the wave strength satisfies 5 < Sq, then in any time interval [h,T] with < h < T , there 
exists a positive constant Sq = eo{S,h,T), such that if the Knudsen number e < Eq, then the Boltzmann 
equation admits a family of smooth solutions f^'^{t,x,£,) satisfying 

sup ||r'''(t,a;,0-M,^£^Q,(t,x,OllL?(^) < Ch,T £*|lne|, 

where T,h,T — {{iix)\h < t < T,\x\ > h,\x — S'it\ > h}, the norm \\ ■ ||j2^ ^^ ^ is \\ ||l|(R3) '^'^'^ 
positive constant Ch,T depends on h and T but is independent of e. Consequently, when e — > 0+ and then 
h — >■ 0+, T — >■ +00, we have 

|ir^''(e)-M[^^^g](OIL.(_^)(i,.T)^0, a.e. in R+ x R. 

Remark 1 Theorem \l.l\ shows that away from the initial time t = Q, the contact discontinuity at x ~ 
and the shock discontinuity at x — s^t, for small total wave strength 5 < 5q and Knudsen number 
e < £o, there exists a family of smooth solutions f^'^(t,x,£^) of the Boltzmann equation which tends to the 
Maxwellian Mjy ^ gj (t, a;, ^) with {V ,U ,Q){t, x) being the Riemann solution to the Euler equations as a 
superposition of a 1-rarefaction wave, a 2-contact discontinuity and a 3-shock wave when e — > with a 
convergence rate e s | In e| . Note that this superposition of waves is the most generic case for the Riemann 
problem. Similar results hold for any other superpositions of waves by using the same analysis. 

Remark 2 The proof of the above theorem crucially depends on the introduction of two kinds of hyperbolic 
waves. The hyperbolic wave I was constructed by Huang- Wang- Yang \21}/ for the compressible Navier- 
Stokes equations to recover the viscous terms to the inviscid approximation of rarefaction wave pattern 
where the rarefaction wave structure plays an important role in the construction. 

The hyperbolic wave II is constructed to remove the error terms due to the viscous contact wave 
approximation. Note that the construction of the hyperbolic wave II can not be done simply around 
the contact wave approximation as the hyperbolic wave I for the rarefaction wave. Otherwise, the wave 
interaction terms thus induced will lead to insufficiently decay in term of the Kundsen number. Instead, 
it is constructed around the superposition of the approximate 1-rarefaction wave, the hyperbolic wave I, 
the 2-viscous contact wave and the 3-shock profile as a whole. Moreover, it also takes care of the non- 
conservative terms in the previous reduced system so that energy estimates can be taken for anti- derivative 
of the perturbation. 

Remark 3 Note that the analysis can also be applied to the vanishing viscosity limit of the one dimen- 
sional compressible Navier-Stokes equations. In fact, the vanishing viscosity limit of the one dimensional 
compressible Navier-Stokes equations in some sense can be viewed as a special case of hydrodynamic limit 
of Boltzmann equation to the Euler equations by neglecting the microscopic effect. 

Remark 4 // the total wave strength S |(f+ — m+ — m_, 0+ — 0-)\ < Sq, then from the wave curves 
defined m (fOS]) . (fr23l) and (fTM]) . we know that (5^S 5"^^ , < cSa where S^\S'='^ ,6^^^ are the wave 
strengths of rarefaction wave, contact discontinuity and shock wave, respectively. 

Let us now review some previous works on the hydrodynamic limits to the Boltzmann equation. For 
the case when the Euler equations have smooth solutions, the vanishing Knudsen number limit of the 
Boltzmann equation has been studied even in the case with an initial layer, cf. Caflisch j6], Lachowicz 
[26] . Nishida ^3] and Ukai-Asona [37] etc. However, as well-known, solutions of the Euler equations in 
general develop singularities, such as shock waves and contact discontinuities. Therefore, how to verify 
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the hydrodynamic limit from the Boltzmann equation to the Euler equations with basic wave patterns 
becomes a natural problem in the process to the general setting. In this direction, with slab symmetry, 
as mentioned earlier, there were studies on each individual wave pattern. For superposition of different 
types of waves, to our knowledge, there is only one result given in |20j about the superposition of two 
rarefaction waves and one contact discontinuity. 

On the other hand, for the incompressible equations, there are works, such as those by Bardos-Golse- 
Levermore, Bardos-Levermore-Ukai-Yang, Bardos-Ukai, Golse-Saint Raymond, Levermore-Masmoudi and 
Sone which studied direct derivations of the incompressible Navier-Stokes equations in the long time scal- 
ing, about which more is known, cf. [21 IH [31 [HI [251 133 [36] and the references therein. In particular, 
Golse and Saint-Raymond showed that the limits of suitably rescaled sequences of the DiPerna-Lions 
renormalized solutions to the Boltzmann equation are the Leray solutions to the incompressible Navier- 
Stokes equations. However, even in this aspect, the uniqueness and regularity of the solution are still big 
issues. Since we will concentrate on the compressible Euler limit in this paper, we will not go into details 
about the incompressible limits. 

Furthermore, the Boltzmann equation provides more information than the classical fluid dynamical 
systems so that it describes some phenomena which can not be modeled by using the classical systems, 
such as Euler and Navier-Stokes equations. This kind of interesting phenomena, such as the thermal 
creep flow in a rarefied gas was known since the time of Maxwell. Some mathematical formulations 
and numerical computations on the basis of kinetic equations were studied since 1960s, cf. the works 
by Sone [35[ 136] . However, the justification of this kind of fiuid dynamics is almost open with rigorous 
mathematical theory. 

Finally, we briefly outline the proof of the theorem. Firstly, we define the individual wave profile. 
Then we introduce the first family of hyperbolic wave by linearizing around the approximate rarefaction 
profile and by adding the viscosity and heat conductivity terms induced by the profile. Then we define 
the first approximation of the superposition of this hyperbolic wave together with the three basic wave 
patterns so that it takes care of the hyperbolicity of the rarefaction wave in the viscous setting. 

Based on this, we linearize the fluid system around this profile and consider the propagation of the 
extra error due to the contact discontinuity approximation and then define a second set of hyperbolic wave. 
By adding these two sets of hyperbolic waves to the superposition of the three basic wave profiles, we will 
perform the energy estimate on the Boltzmann equation with suitable initial data through the macro- 
micro decomposition. Precisely, for the macroscopic component, we will consider the anti-derivative of 
the perturbation after applying a hyperbolic scaling. By using the dissipation in the fluid- type system 
and the linearized Boltzmann operator on the microscopic component, we can close the energy estimate 
through a suitable chosen a priori assumption. Then the statements in the theorem follow. 

The rest of the paper will be arranged as follows. In Section 2, we will construct the approximate 
solutions to the Boltzmann equation corresponding to the basic wave patterns to the Euler system. Then 
we obtain the detailed information on the difference between the Riemann solution to Euler system and 
the approximate solution to the Boltzmann equation by the construction. In Section 3, we will construct 
a family of solutions to the Boltzmann equation around the approximate solution by using energy method 
to close the a priori estimate. Since the proofs of two Propositions 13. II and 13.21 about the lower and higher 
order energy estimates respectively are very technical and long, we put them to the Appendices. 

Notations: Throughout this paper, the positive generic constants which are independent of e, T, h 
are denoted by c, C, Ci(i = 1, 2, 3, • • • ), while Ch,T represents a generic positive constant depending on 
h and T but independent of e. And we will use || • || to denote the standard L2(R]dy) norm, and 
II • ll^i (i = 1, 2, 3, • • • ) to denote the Sobolev i?*(R; dy) norm. Sometimes, we also use 0(1) to denote a 
uniform bounded constant which is independent of £, T, h. 

2 Approximate Wave Patterns 

In this section, we will construct the approximate wave profile that consists of three basic wave patterns 
and two hyperbolic waves. For this, we will firstly recall the construction of the approximate rarefaction 
wave for the Boltzmann equation. Then we will introduce the hyperbolic waves I to correct the error 
terms coming from the hyperbolic approximation. Then we will construct the viscous contact wave 
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to Boltzmann equation and study the non-conservative error terms. The viscous shock profile to the 
Boltzmann equation wih then be recalled. With the above wave patterns, we will introduce the hyperbolic 
wave II to take care of the error terms due to the viscous contact wave by avoiding the interaction between 
the viscous contact wave with the wave patterns defined earlier. 



2.1 Rarefaction Wave 



For the rarefaction wave, since there is no exact rarefaction wave profile for either the Navier-Stokes 
equations or the Boltzmann equation, the following approximate rarefaction wave profile satisfying the 
Euler equations was introduced in |32[ 139] . For the completeness of the presentation, we include its 
definition and the properties obtained in the above two papers as follows. 

If {v-,ui-,9-) € then there exists a 1-rarefaction wave {v^'^ ,Ui^ , E'^'^){x/t) which is 

a global solution to the following Riemann problem 



Vt - Uix = 0, 
Wit +Px = ^, 

Et + (pui)^ = 0, 



{v+,ui+,e+),x > 0. 



(2.1) 



Consider the following inviscid Burgers equation with Riemann data 

Wt + WWx = 0, 



w(t = 0, x) 



w+, 



X <0, 
x>0. 



(2.2) 



If w^ < wjf- , then the above Riemann problem admits a rarefaction wave solution 

W-, J < w^, 
f, w-<f<w+, 



,x , 



w'-{t,x) = w^{-) 



(2.3) 



As in [39| , the approximate rarefaction wave (y^i,[/^i,e^i)(<,a;) to the problem can be con- 
structed by the solution of the Burgers equation 



Wt + WWx = 0, 

w(Q,x) = Wa{x) = w{-) = — 



W+ — W- x 

tanh — , 



(2.4) 



where cr > is a small parameter to be determined later to be £ 5 . Note that the solution wJJ.(t, x) of the 
problem (|2.4|) is given by 

wl{t,x) = Waixo(t,x)), X = XQ{t,x) + WaiXo{t,x))t. 

The smooth approximate rarefaction wave profile denoted by {V^'^ , U^^ , &^^){t, x) can be defined by 



' S^''{t,x) ^ s{V^'Ht,x),e'^'{t,x)) 
w± ^ Ai± := Ai(w±,6'±), 
w:{t,x)^X,{V''Ht,x),s+), 

U^'{t,x)^ui+ 



(2.5) 



Xi{v,s+)dv, 



Ui\t,x) = 0, i 2,3. 
Note that {V'^\U^-\Q^'^ ){t, x) defined above satisfies 



Vt""' ~ u?- - 0, 
t/it + P^' = 0, 
U^'=Q, * = 2,3, 



(2.6) 
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where P^^ = p{V^^ ,<d^^) = I^ttt and = + ■2'^ ■ The properties of the rarefaction wave 
profile can be summarized as follows. 

Lemma 2.1 (i39^ ) The approximate rarefaction waves (y'^\U^\Q^^){t,x) constructed in (^3)) have 
the following properties: 

(1) U";^ (t, a;) > for x G R, t > 0; 

(2) For any 1 < p < +00, the following estimates holds, 

|!(F«Sf/fSe«0.|Und=.) <t^niin{(5^ia-i+i/P, (5«i)i/Pi-i+i/P}, 
||(l/«SC/i«\e«Ox.||Lp(d.) <C^min{5«ia-2+i/P, ^-i+i/p^-i}^ 

where the positive constant C depends only on p and the wave strength; 

(3) Ifx> X^'p, then 



\dUV''\U^\e^^){t,x)\ < -^e-^^, fc = l,2; 
(4) There exist positive constants C and ctq such that for a G (0, (To) '"^'^ t > 0, 

sup I (y^i , , ) (i, x) - (w'-i , u'-i , (-) I < - [fT ln(l + i) + ct| In a|] . 

2.2 Hyperbolic Wave I 

Since the whole wave profile consisits of a shock wave whose rate of change in the shock region is of 
the order of we have to consider the anti-derivative of the perturbation in order to cope with the 
correct sign as in the stability analysis. From p.6p . we know that the approximate rarefaction wave 
x) satisfies the compressible Euler equations exactly without viscous terms. Thus if 
we carry out the energy estimates to the anti-derivative variables, the error terms due to the viscous 
terms from the approximate rarefaction wave are not good enough to get the desired estimates. In order 
to overcome this difficulty, we introduce the hyperbolic wave I to recover these viscous terms. 

This hyperbolic wave denoted by (di, (i2, d3)(i, x) can be defined as follows. Consider a linear system 

d\t - dij: = 0, 

d2t + {P^^d, +p^ld, +pl^d,y, = -e C yl (2.7) 



d^t + [{pui)^'di + {pu^)^ld2 + {pui)fd:i\^ - e{ " + -e{^^ i^). 



where p = ^ = p{v,u,E) = and p^^ = py{V^\U^\£^^) etc. Note that the left hand side of 

the above system is the linearization of the Euler equation around the rarefaction wave approximation. 
We want to solve this linear hyperbolic system (|2.7p on the time interval [h,T]. Firstly, we diagonalize 
the above system by rewriting it as 




t 



A 





(2.8) 



where iff ^ = .(Me;^^)^, ^ .(..(e^^e^^ Here, the matrix 

0-10 

A^^^i P^^ pI^ 

{pu,)^^ {pm)^l {pm)l^ 
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has three distinct eigenvalues Af^ := Xi{V'^\s±) < = Af^ < X3{V^^,s±) := Af^ and the correspond- 



ing left and right eigenvectors denoted l^^ ,r^^ {j = 1,2,3) respectively, satisfy 

LRiA^^r"^ ^ diag(A^,0,A^^) ee A^S L^'R"' Id., 

Here L^^ = (/f ^ , l^' , l^^' )*, = (rf ^ , ^ , rf ^ ) with /f ^ = , J/f ^ , s+) and rf ^ = r,(T/^i , C/f ^ , s+) (i 

1,2,3) and Id. is the 3x3 identity matrix. Now we set 



{Di,D2,D:iY ^L"^{di,d2,dz) 



Then 



(di,d2,d3)* =-R'^^(i^l,i^2,^3)*, 

and Z?2, -D3) satisfies the system 



(2.9) 
(2.10) 



\ /A 

(2.11) 

Due to the fact that the 1— Riemann invariant is constant along the approximate rarefaction wave curve, 
we have 




Li 



Ri 



-AfLf^ 



Substituting the above equation into (|2.11l) . we obtain the diagonalized system 



D2t = b^^H^' + b^^H^' + a^^V^'D2 + a^^V^'D^, 



(2.12) 



where af^ , bf^^ are some given functions of V^^ , U^^ and 5*^^ = s± . Note that in the diagonalized system 
(|2.12p . the equations of -D2, are decoupled from Di because of the property of the rarefaction wave. 

Now we impose the following boundary condition to the above linear hyperbolic system (j2.12p in the 
domain (t, x) e [h, T] x R: 



Di{t = h,x) = 0, D2{t = T,x)^ Dsit = T, a;) = 0. 



(2.13) 



With this boundary condition, we can solve the linear diagonalized hyperbolic system (|2.12p under the 
conditions (|2.13p . Moreover, we have the following estimates on the solution. 



Lemma 2.2 There exists a positive constant Ct.T independent of e such that 
(1) 

II Olli^Cd,) <C;,,T i = 1,2,3, fc = 0,1,2,3. 

(2) If X > Xi+t, then we have 



1 l^-Ai+t 

\d,{x,t)\ < Ch,T -e ^ 

(7 



\dix{x,t)\ < Ch,T -^e" 



-, z = l,2,3. 



The proof of Lemma [2.21 can be done similarly as in [21] for the compressible Navier-Stokes equations. 
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2.3 Viscous Contact Wave 

In this subsection, we construct the contact wave {V^^,U^^,e'^^){t,x) for the Boltzmann equation 
motivated by [351. Consider the Euler system (|1.16p with a Riemann initial data 

iv,u,e)it^O,x)^{ (2.14) 

where u± — (ui±,0,0) and v± > 0,9± > 0,ui± are given constants. It is known (cf. [34j) that the 
Riemann problem (|1.16p . (|2.14p admits a contact discontinuity 

«u^^r'^)(t,.) = | (--'--'^-)' ^<0, (2.15) 
^ ' ' ^ \ {v+,u+,0+), x>0, ^ ^ 

provided that G CD{v+,ui+,9+), that is, 

S^*- 26'+ 
= P- ■= o — =P+-=o — ■ (2-16) 



Then for the Navier-Stokes equations, by the energy equation (ll.2ip zi and the mass equation (|1.2ip i with 

3p+ ■ 



cf. [20], we can obtain the following nonlinear diffusion eciuation 



% = e{a{9)e,U a{e) = ^P±^. (2.17) 



From [T and |llj . we know that the nonlinear diffusion equation p.l7p admits a unique self-similar 
solution 8(77), 7] = , ^ satisfy 

then 8(t, a;) has the property that 



solution 8(77), 7] = ^ {i+t) ^^tisfying the boundary conditions 8(±oo.t) — 9±. Let (J*^"^ 



e.(^,x)^ / e-^, (2.18) 

with some positive constant c depending only on 9± . 

Correspondingly, we can define the Navier-Stokes profile by 

^^f8, 

. 2sa(8). . (2.19) 

3p+ 

For the Boltzmann equation, if we still use the above Navier-Stokes profile {V, U, 8), we can not get any 
decay with respect to the Knudsen number e due to the non-fluid component. Hence, we construct a 
Boltzmann contact wave as follows. Set 

= I^^M {PiKi(^9^0?^ +C • U^M]}, (2.20) 

and 

nfi^ = L^[ei~^G^^ + iPi(aG.^^)) - Q(G^^,G^^)], (2.21) 

where (V^^ , U'~^^ , <d'~^^){t, x) is the viscous contact wave for the Boltzmann equation to be constructed 
later. 

Note that for the Boltzmann equation, the leading terms in the energy equation ([TT21J4 can be written 

as 

9, +p+m. = si^y. - 1 ieil^pnf,^de + ^i+/^e?nf,^de (2.22) 

By the definition of Hf^^ in ([2^, we have 

- / ^aiePnfi^de + u,+ J i^^nfi^dC = An + Ai2, (2.23) 
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where 



All 



with gii — gii{v, u, 9), {i ~ 1, 2, 3, 4) being smooth functions of (w, u, 0), and 

= o{i)e' [{\v,\ + \u,\ + \e,\ + le^^i + |c/f ^|)|c/f ^1 + \u,\\e^^\ + \ug^\ 

Thus, by choosing the leading term and dropping the higher order term Ai2x in (I2.22p . we have 



■.)x + -r^llx 

5 



(2.24) 



(2.25) 



(2.26) 



where a{9) is defined in (|2.17p and An is defined in (I2.24p . To represent the microscopic effect on the 



wave profile, we want to define Q'-'^ to be close to 9(- 



= ) + Q"^{t,x) with 9 being determined 



by (|2.17p . (|2.18p and represents the part of the nonlinear diffusion wave coming from the non-fluid 
component not appearing in the Navier-Stokes level. Moreover, the term Q"'^ decays faster than O so 
that it can be viewed as the perturbation around the Navier-Stokes profile Q. To construct Q"-^ , we 
linearize the equation (|2.26p around the Navier-Stokes profile Q and drop all the higher order terms. 
This leads to a linear diffusion equation for S>"-^ 



ft"/ 



e 



(a(e)esO. + £(«'(0)0xe"-^). 



rAn 



(2.27) 



where An = £^(511 + ^gi2 + gi?){^xf + e^gu^xx with gu = gu{V , U, 6) {i = 1, 2, 3, 4). Integrating 



(|2.27p with respect to x yields that 



where 



^i(i,x) 



Note that An takes the form of -^A"^ 



e"^(i,x)dx. 



) and satisfies that 



(2.28) 
(2.29) 



|An| = 0(5'^-°)e(l +<)"'e ^"Ci)-^**) , as a; -> ±00. 



We can check that there exists a self-similar solution ^i(- 



= ) for p.27p with the boundary conditions 



E!i(— 00) = 0, Si(-l-oo) — Here Si4_ can be any given constant satisfying < 6'-^^ . It is worthy to 

point out that even though the function Si(i, x) depends on the constant 8"/(t, a;) = 'B.ix{t, a;) — > 
as a; — ^ ±00. That is, the choice of the constant Si_|_ has no influence on the ansantz as long as |Si+| < 
^CD ^ From now on, we fix Si+ so that the function Si(t,a;) is uniquely determined and its derivative 
Sia; = 6"/ has the property 



|e"/| = |Si,| = 0(<5^^)e5(l + <)-i 



4a(e±)e(l + t) 



as a: — )■ ±00. 



(2.30) 



Then we apply the similar procedure to construct the second and the third components of the velocity 
of the contact wave denoted by U^^ {i = 2, 3) as follows. The leading part of the equation for Ui in 
{i = 2,3) is 



e — — — U,: 



Firstly, we have 



where 



A.i = 



29 



I' 



(2.31) 

(2.32) 
(2.33) 
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with gij, {i = 2, 3, j = 1, 2, 3, 4) being the smooth functions of (u, u, 9) and 

Thus we expect the viscous contact wave Uf^ {i — 2, 3) to satisfy the following linear equation 



"-'it 



CD _ ttCD\ 



26 



Aiia;, i = 2,3, 



(2.34) 



(2.35) 



where where = £^(3^1 + 3|^gj2 +5»3)(0a;)^ + e^5j46xx with g^^ = gij{V,U,Q) {i = 2,3, j = 1,2,3,4). 
Integrating (12. 35^ with respect to x yields that 



26 



where 



Note that A^i takes the form ■j^A'^{ 



(2.36) 
(2.37) 



— 2,3 and satisfies that 



I A,i I = 0{6^°)e{l + ty^e' 4a(«±%(i+t, ^ as x ^ ±c3o. 



We can check that there exists a self-similar solution 



=) for (|2.35p with the boundary con- 



ditions Si(— 00) = 0,Si(-|-oo) = {i = 2,3). Again, here Si+ can be any given constant satisfying 

jSi+l < S"-^^ . As we explained before, the choice of the constant Si+ has no influence on the ansantz as 
long as < (5*^^. We fix Si+ so that the function Si(i, x) is uniquely determined and the derivative 

Six = has the property 



CD I 



0((5^^)£5(i + i)-2e 45(.±).(i+t) ^ ^ x^±oo, 



(2.38) 



with 6(0±) = max{a(0±), ^^^f^)}. 

In summary, the viscous contact wave [V'-^^ , U'-^^ , Q'~^^){t,x) can be defined by 

c/f = ui+ + — U(e)e, + e(a(e)e"Ox + ^An 
3p+ L 5 

C/P^ = S„, (z = 2,3), 

e'^-^ = e + 9"-^ + H, 



(2.39) 



where 



H=- 



eV 



CD 



P+ 



a(e)et+(a(e)e"0 



46^(6) 
3p+ 



(a(e)e,), + (a(e)e"/)..J-2^'''' Uf^nd^^ (2.40) 



is chosen such that the momentum equation that the viscous contact wave satisfies has an error term 
with sufficient decay in e. Without H, the error term in the momentum equation decay like £2 . In order 
to get e order decay in the error term, we should introduce the higher order approximate term H in the 

definition of in the contact wave. Here, J ^^Ilnd^ in (|2.40p is the corresponding function defined 
in ((2?2T]) by replacing both the variables {v, u, 9) and {V^^, U^^, Q^^) by {V, U, 9), and it satisfies 



/ 



C?nnde = 0(1)6^1(92, e,,)|. 



(2.41) 
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Hence, by (I^JS)) . (j^:^ and (I^TiT]) . we have 

i? = 0(<5'^'°)e(l + t)~2e"^fi^, as a; ^ ±cx). 
Now the contact wave {V^^, C/*^^, Q'^^){t,x) defined m (|2.39p satisfies the following system 

rCD jjCD 



(2.42) 



r ^ - C/rT^ 0, 



CD , pCD ^ f£/ 



3 ' VCD 

yCD 
CDtjCD 



3 



yCD 

2^^^ yCD I'-' 



yCD 



(2.43) 



where P^^ = f^, = 9^^ + - j Ci^n^df {^ = 1,2,3) and - j 6^nf,^d^ are 

the corresponding functions defined in p.2ip by replacing the variables {v,u,6) by (^'^■°, C/*^^, 0*^^), 
respectively. Moreover, 



^1 "5^^"* "15^^ ^-^yV C/i. j^ + y ei(nn -IlnU^ 



2e 



as a; — > ±oo, 



(2.44) 



and 



Q 



CD 



/ Me^^)~A^(e) 



as X ±oo, i = 2, 3, 



5e 

y 



(a(e^^) - a(e) - a'(e)(e^^ - e)) (e, + e^O + a(e^^)i?, 
3p+e f K{ec^)&(^^H 



2Ug°H 4e^(e^^ 



3V^CD 



3Ft^^ 



2 Ve<^-D(0CD_^) 

i=2 



i=l 

as X — >■ ±cx). 



(2.45) 



(2.46) 



with some positive constant c > depending only on 9± and Ai2x (i ~ 1,2,3) being the corresponding 
functions defined in and ([234)1 by replacing both (w, u, 6*) and (V^^^, C/<^^, O*^-") by (f , {/, 9). 

Note that from (|2.18p . we have 



|(y^^,i7^^,e^^)(t,x) - {v''',u^^,e^'')it,x)\ = 0(l),5^^e~5?Ti+i). 



(2.47) 



2.4 Shock Profile 

In this subsection, we will firstly recall the shock profile F^^{x — sst,^) of the Boltzmann equation (|l.ip 
in Eulerian coordinates with its existence and properties given in the papers by Caflisch-Nicolaenko [7] 
and Liu-Yu |30j , |31) . And then we will state the corresponding properties in the Lagrangian coordinates 
used in this paper. 
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First of all, F^'^^ (x — s^t, ^) satisfies 



-S3 



(2.48) 



where ' = — x — s^t, u± = (ui±,0,0) and {p±,u±,9±) satisfy Rankinc-Hugoniot condition 

-S3{p+ - P-) + {p+ui+ - P-U1-) = 0, 
-S3(p+Mi+ - p-ui-) + {p+uj_^_+p+ - p-ul_ = 0, 

-Sz{p+E+ - p-E_) + (p+ui+£'+ + P+UI+ - p_ui_£'_ - p-ui-) = 0, 



and Lax entropy condition 



(2.49) 



(2.50) 



with S3 being 3-shock wave speed and Af" = ui + -^^-p^ being the third characteristic eigenvalue of the 

Euler equations in the Eulerian coordinate and A|lj_ = ui± + 

By the macro- micro decomposition around the local Maxwellian M'^^ , set 



where 



with 



m^'-{x,t,C) = M 



F^^ {x, t, = (x, t, + G^^ {x, t, , 

p^^{x,t) 



Ip^3 ^u^a fiSs-^ 



{x,t,0 



y^{2TTR9S3{x,t))^ 



pSa 

p^'uf' 



I = / I F^'^i^^t^Odt « = 1,2,3. 



(2.51) 



2 ' / \ 2 

With respect to the inner product (•,-)mS3 defined in (|1.6p . we can now define the macroscopic 

s • • • ■ s 

projection Pq^ and microscopic projection Pj^'^ by 



Po'9 = )M-3xf , Pf^9 = 9 - Po^ff, 

i=0 



(2.52) 



where x*^ (j = 0, 1, 2, 3, 4) are the corresponding pairwise orthogonal base defined in (|1.7p by replacing 
(p, w, 61, M) by (p^^' , u^-' , 6*^3 , M^^ ) . 

Under the above macro-micro decomposition, the solution F^^ = F^^{x — sst,^^) satisfies 

and the Boltzmann equation (|2.48p becomes 

(M^-^ +G^^)i +ei(M^-^ +G^-^)^ = i[2Q(M'5-\G'5^') + 0(G'5-%G^^')]. (2.53) 
Correspondingly, we have the following fluid-type system for the fiuid components of shock profile: 
f pf^ + (p^^iif^). = 0, 

{p'^uf^t + [p'Hu'!'r+p'% = ^{pie'^ut)^ - J f?nf3d5, 
{p'^uf^)t + {p^^uf^uf% = eipie'^ui^)^ - J CiC.nf^d^, z = 2,3, 

[p^3(6''^3 -t- 



3 



(2.54) 
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In fact, from the invariance of the equation (|2.48l) by changing ^ with — ^ and the fact that Ui± = 0, we 
have uf' = J fi^nf^df = for i = 2, 3. 

And the equation for the non-fluid component G'^^ is 

Here Ljyjsa is the hnearized colhsion operator of Q{F^^ , F^") with respect to the local Maxwellian M'^=': 

L^s.g = 2Q(M^^5) = Q(M^^5) + Qig,M'')- 

Thus 

G^^ =eLj^i.3[Pf (ClMf^)]+^^^ 



Now we recall the properties of the shock profile F^^ {x — s^t, ^) that are given or can be induced by 
Liu-Yu in Theorem 6.8, |31) . 

Lemma 2.3 (131] ] If the shock wave strength 5^^ is small enough, then the Boltzmann equation (jl.ip 
admits a 3-shock profile solution F^^ {x — s^t, f ) uniquely up to a shift satisfying the following properties: 

(1) The shock profile converges to its far fields exponentially fast with an exponent proportional to the 
magnitude of the shock wave strength, that is 



-P±,wf^ -'"i±,6'^' -6'±)l < C'(5^'e-'=±— as i? ^ ±cx), 

as 1? — >■ ±00, 



with 5^^ being the 3-shock strength and Mo being the global Maxwellian which is close to the shock 
profile with its precise definition given in Theorem 6.8, \31f . 



(2) Compressibility of 3-shock profile: 



(Af)^ < 0, Af - + 



(3) The following properties hold: 



„S3 „,S3 aSs f\E\ 



Mn 



where A '-^ B denotes the equivalence of the quantities A and B, and 
uf'=0, J a^n^^d? = 0, z = 2,3, 

(xS3\k-l 



KI^I)|ggG^3| 
Mo 



Mo 



-d^)\ k>l, 



I J ^lV^mfd^\<C6''^\uf'^\, 1 = 1,2,3,4, 
where ipi{^) {i = 1,2,3,4) are the collision invariants defined in (|1.4I) . 
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Now we rewrite this shock profile in Lagrangian coordinate by using the transformation (jl.l4l) and use 
{t,x) for the Lagrangian coordinate to distingish it from the Eulerian coordinate (t,x) at this moment. 
Then the shock profile in Lagrangian coordinate can be written as F^^ {x — s^i, ^) with S3 determined 
by the 3-shock wave curve given in (|1.24p . First, from the Rankine-Hugoniot condition in Eulerian and 
Lagrangian coordinates, we have 

-S3{P+ - P--) + {p+ui+ - P-Ui-) = 0, 

and 

-33(1;+ - V-) - - Ml-) = 0, 

with v± — respectively. Thus we have the following relation between S3 and S3 

S3 = P±(s3 - Mi±)- (2.55) 
On the other hand, we have from (|2.54P -^ that 

p^''{x - S3t)[s3 - uf^(2; - S3t)] = const. = p±(s3 - ui±). (2.56) 



Note that 



ft pX pt 

x-ssi = p^^uf''{-S3T)dT + p^" {y - S3T)dy - P±{s3 - ui±)dT 
Jo, Jo__ Jo^ 

p^'^(y-S3T)dy- f p^'{-S3T)[s3~uf'{-S3T)]dT 



- 0_ 

t px—s^t 



- / p'^'uf'{^7S3T)dT + 

Jo Jo^ 

= P^'{y-S3T)dy-S3 p^^{~S3T)dT = I p^'iv)dv 
Jo Jo Jo 

where in the second equality we have used (|2.56l) . 

This shows that under the Lagrangian transformation (|1.14p , the shock profile F^'-^ {x — s^t, ^) in Eule- 
rian coordinate can be exactly transformed to the shock profile F^^ (i — 53^, f ) in Lagrangian coordinate. 
Moreover, we have 

F^HV,0=P''F^'{^,0 

with rj = x — s^t. 

For simplicity of the notations, from now on, we use {t, x) to denote the Lagrangian coordinate and 
F^^ {-q, £_) with •q = X — s^t to denote the 3-shock profile of Boltzmann equation in Lagrangian coordinate. 
And in the Lagrangian coordinate, we have the following Lemma. 

Lemma 2.4 Assume that {v-,u^, 6^) e S3{v+,u+,9+), then there exists a unique shock profile F^^{ri,£^) 
with Tj ^ X — s^t up to a shift, to the Boltzmann equation (ll.lSp in Lagrangian coordinate. Moreover, 
there are positive constants c± and C such that for 77 e R, 

S3Kf-' = -C/f,^ > 0, 
C/f^ =0, j 66nf de = 0, i = 2,3, 

{\V^'' - v±\,\U^'' - ui±\,\Q^^ - 6±\) < CT^^e- '^^""' , as 77 ^ ±00, 



Furthermore^ we have 

Mo 

and 

|a,^(F^^f/f^e^^)l < l(Kf^t/f,^e,t)l' ^ > 2, 



Mo - ^ J Mo 
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and 



S3 IrrSa 1 



with ipiiS,) being the collision invariants. 



* = 1,2,3,4, 



Furthermore, we have 



3^ 



(2.57) 



1=2 



yS3 



where £^ = <d^ + - — j-'- and {v± , m±, 0±) satisfy Rankine-Hugoniot condition and Lax entropy condition 
and S3 is 3-shock wave speed. 

Correspondingly, we have the following equation for the non-fluid part of 3-shock profile. 

Gf^ - ^Gf-' + A.pf3(^iMf3) + _Lpf3(^^Gf^) - i [LmS3G^-^ -t- Q(G^3^ G^^')] . 

Here, Lj^sg is the linearized collision operator of Q{F^^ , F^^) with respect to the local Maxwellian M'^^: 

LMS3.g = 2Q(M^%g) = Q(M^%<?) + Q{9,M'')- 

Thus 



(2.58) 



S3^ 



2.5 Hyperbolic Wave II 

The purpose of this subsection is to construct the second hyperbolic wave. Up to now, we can define the 
following approximate composite wave profile {V,U,£){t,x) 



V 

Ui 
£ 



{t,x) 







Vv* \ 












-E* 1 



(2.59) 



[7, = !7f^,i = 2,3. 



where £ = Q + {V^^, U^^ ,£^'^){t,x) is the 1-rarefaction wave defined in (|2.5p with the right state 
{v+,uij^,E^) replaced by (1;*, £'*), {V^^ ,Ui^ ,£'~^^){t,x) is the viscous contact wave defined in 
(j2.39p with the states {v^,ui^,E^) and (w+,mi+,£'+) replaced by (w* , wi*, E'*) and {v* ,u\,E*) respec- 
tively, and {V^^ , [/f ^ ^£^^){t,x) is the fiuid part of 3-shock profile of Boltzmann equation defined in (|2.57|) 
with the left state {v-,ui-,E-) replaced by {v* ,u\,E*). 
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Moreover, we can check that this profile satisfies 



Vt - C/i. = 0, 

_ 4 

,^i{Q)Ulx. 



e.^'iFxdi 



V 



U^^lFxd^ 



CD 



2,3, 



8t + (PC/ija; = e( — — )x + 77 ):e + Z^^( f> 



V 



V 



V 



where P = 8), Qf^ (i = 1,2,3,4) are defined in (P^ii)) . (piiSl) and (pTig)) . respectively, and 



P - {p^'di + p^ld2 + p^E'd^) ~ - P^^ - P'^^ +p* 

3 L y yJ^i yc^" 



■V V 



Q4 



k{q)Qx K(e«i)0^' K(e^^)e^^ K(e^3)ef3- 



V yc^^ F'^^ 

3 L y yfli yen 

^ rM(9) /i(e^^) 



j=2 



From (|2.60p . we have that 



e= +d3 + e^-° + e^^ - (61, + r) 

_1 [|C7|2 „ |C;«.|2 _ pCD^2 |f^53|2 ^ ^2^ ^ (^*)2j^ 



Direct computation yields 



Qi 



CD 



CD 



til* 



e"^" ~ e,,v'^'' ~ V* ,u^' - ui,e^' ~ e*)\ 
+\{v^' - V* , t/f ^ - , e^^ - 0* ) 1 1 (F^^ - V* , c/f ^ - ui* , t/f ^ - , e^^ -e*)\ 

+e\U^^'\\{V^'' -v,,e^'' -0,,V'^' -v^Q"^' -0*)\ 

+e I Ug'' 1 1 (F^i - i;. , - , di , d2 , d3 , T/^^' - i;* , O'^^ - r ) I 

+ 0{l)[\{dud2,d3)\^ +e\d2x\+e\U{l'\\{dud2,d3)\ 

■~ Qii + Qi2i 
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and 



+£ I (L/g , ef ^ ) 1 1 - z;,, e^^ - e„ , d2 , 4 , ~ v* , c/f ^ - , e^^ - r ) | 

+0(l)[|(di,d2,d3)l'+£|(d2.,d3x)l+e|(C/i'^\e^^)||(di,d2,d3)| 



(2.64) 



Here, Qn, Qi (i = 2,3) and (54i represent the interaction of waves in different families, Q12 and Q42 
represent the error terms coming from the approximate rarefaction wave and the hyperbohc wave I. 

Firstly, we estimate the interaction terms Qn, Qi {i = 2,3) and Q41 by dividing the whole domain 
il = {{t,x)\{t,x) G [h,T] X R} into three regions: 

={{t,x) en\2x<XiJ}, 

^CD ^ {(t,x) e^\Xi*t <2x < sst}, 

ns, ^{{t,x)en\2x>s3t}, 

where Ai* = Xi{v^,9^) and S3 is the 3-shock speed. 

From Lemma |2.1[ we have the following estimates in each region: 



• In rtn^ , 



CD 



CD 



• In ilcn, 

|(l/«i -«*,di,d2,d3)| =0(l)e- 



2lan-Ai.t| 



0(l)e- 



|y^^3_„*| = 0(l)e-^^%^ = 0(l)e- 



• In 17 



S3' 

|(l/«i -t;*,di,d2,d3)| =0(l)e 



0{l)e- 



\{V'^^ -v^,V'^^ -v*)\ = Cl(l)(5C^e"'8^i+*)' e"SriTTT = 0(l)e- 



Note that we just give the pointwise estimates of ^component and di {i = 1, 2, 3) in each region, similar 
estimates hold also for the Ui and Q components. In summary, we have 



KQii, Q2, Q3, Q4i)| = Ch,T e 



(2.65) 



with a ^ and for some positive constants Ch.T and Ch independent of e. 

In order to remove the non-conservative error terms Qf^, {i = 1,2,3,4) coming from the definition 
of the viscous contact wave, we now introduce the following hyperbolic wave 6 = (61, 621, ^22, ^23, ^3): 



bit - hix = 0, 

&21t + [Pvbl + PuM + PuM2 + ^'n3^23 + ^'iJ&s]:^ = "Qf^, 
&22t — —Q2^, 
&23t = —Qs^, 

ht + [{PUi^hi + {PUi)u,h2i + {PUi)u,b22 + {PUi)u,b23 + iPUi)Eb3U = -Q 



(2.66) 



CD 
4 I 
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where P = ^ = PiV: U, £) and = Pv{V ^ 0, £), etc. For later use, we denote 62 = (621, ^22, &23)- Now 
we want to solve this linear hyperbolic system (|2.66p on the interval [h, T] . Firstly, we diagonalize the 
above system. Rewrite the system (j2.66p as 





\ 


hi 




622 




623 




V b3 


J 



b2i 

b22 
b23 

\b3 J 



( \ 

V -of^ J 



(2.67) 



where the matrix 



A{V,U,£) 






-1 











Pv 


p 


p 


p 


Pe 

































V (PC/i)„ {PUi)u, iPUi)u, {PUi)u, {PUi)e J 



has three distinct eigenvalues Ai := Xi{V,P) < Q — X2 < X3{V,P) := A3, (here A2 being 3-repeated 
eigenvalues) with the corresponding left and right eigenvectors denoted by 



It holds that 



h,l2i,l22, 123, k; ri, r2i, r22, r23, rs. 

LAR = diag(Ai, A2,A2,A2, A3) = A, 
LR = Id. 



Here L = (^i, ^21, ^22, ^23, ^3)*, -R = (»^i, ?'2i, r-22, ?'23, '^3) with L = L{V,U,£) and R = R{V,U,£) {i 
1,2,3) and Id. is the 5x5 identity matrix. Specially, we can choose 



Set 
then 



= (p,-i7i, 0,0,1), [22^(0,0,1,0,0), [23^(0,0,0,1,0). 

B = {Bi, B21, B22, B23, B3Y — L ■ (&1, ^21, ^22, &237 ^3), 
(&1, 621, ^22, &23, ^3)* = R ■ {Bl, B21, B22, B23, B3Y, 



(2.68) 

(2.69) 
(2.70) 



and B satisfies the system 

/ ^1 \ 

-B21 

B22 + A 

^23 

For simplicity, denote 

Q^^ = (0 
So we obtain a diagonalized system 



/ B, 


\ 




( 









( 






( 


\ 


B21 














B21 






B21 




B22 




^ L 








+ LtR 


B22 




+ L^RA 


B22 




B23 








-QT 






B23 






B23 




V B3 


J 


- X 


\ 


-Q^ 


I 




\ B3 


) 




\ B3 


) 



(2.71) 



(2.72) 



Bit + {XiBi)^ = h ■ + {ht + \hx) ■ nB, + lifY^ r2jB2j, 



1=1.3 



B2U = hi ■ Q^^ + ('"21* + ^»'2ix) • r,B^ + hit ■ Y ^2j-S2j ' 



4=1,3 



^22* — h2 ■ Q 
^23t = ^23 



CD 



(2.73) 



^23* — ^23 ■ Q*^^, 



B3t + (A3B3)^ = ^3 • + J2 iht + Xthx) ■ nB, + hfY ^2j^2i- 

i=l,3 j = l 
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Now we impose the following boundary condition to the linear hyperbolic system (|2.73p on the domain 
{t, x) e [h, T] X R: 

(Bi,B2i,B22,B23,B3)(t = = 0. (2.74) 

We can solve the linear diagonalized hyperbolic system (|2.73p under the condition (|2.74p to have the 
following lemma. 

Lemma 2.5 There exists a positive constant So such that if the wave strength 6 < Sq, then there exists a 
positive constant Ch.T which is independent of e, such that 

621, &22, 623, h)(t, + y II V |C/Sl^(^i, &21, 622, &23, h)(t, ■)\\lHd.)dt 

< Ch,T £^''^ fc = 0,l,2,3. (2.75) 

Proof: From the wave curves defined in (fT^ - lfm)) . we know 5'^^ + 5^-^ < CSq. Let TV = ^ and 
S < Sq 1, then we have 

< CO < (-) , (— ) < Co < 00, 

where V is defined in (|2.19p and co, Co are independent of Sq. 

Without loss of generality, we assume that v* = 1 and > 0. li v* ^ 1, then we just replace V"^^ 
V 

and (y^3)±w by (_)±A^ and ( )±^, respectively. 

V* V* 

Firstly, multiplying the equation (|2J3| i by [(V")"" + {V'^'-'y^]Bi, we obtain 

+ (Ai(F,P)i?i).[f-^ + iV'T'']Bi ^ (2.76) 

= [V~^ + (F^3)-^]Bi (h • + ^ (Fit + hh.) ■ nB, + ht ■ J2 ^2J-B2j 

^ i=l,3 j=l 

Now the second and third terms on the left hand side of (|2.76p can be estimated by 



< Ch,TBl + C{% + |f/g|)B? - ]^N\\i{V,P)\V-''-^%Bl 

+ M:^)(^^3)-A^-i|c;g 1^2 + (. . . 
2 S3 



(2.77) 



Then we estimate the right hand side of ()2.76p term by term. On one hand, 

\[V-N (1/S3)-N]5^;- . gCD| < ^^2 ^ C\Q^°\^. (2.78) 

On the other hand, we have 

\IV-N ^ (^^^3)-^'](^-^^ + Xiri,) • nBf\ < Ch,TBl + CV^BI + C\U^^\Bl, (2.79) 



< C;,,t|S|2 + C\U^^\{Bl + |B2l|2 + IB22P + 1^231'), (2. 
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and 

< CMB! + Bl) + C%{Bl + Bl) + C|Af^ - s:,\\ut\{Bl + 

+C/,,T[e"^e-^ + |(di:.,d2:.,d3:.)|](S? + 

< Ch,T{Bl + + CV^{Bl + S^) + C(5|C/i^'|(b2 + B^). (2.81) 
Substituting (|2.77p - (|2.8ip into (|2.76p and choosing iV large enough give 

(1 [f-^ + {v^^r'']Bi] ~ \n\\,{v,p)\v-''-'\v,\bi - \n{i + M:^)(yS3)-A^-i|c;53|^2 

> - C,,r|g''^P - C\Ut\iz l^2.f - CiVilBl - C5\U^^\Bl 

(2.82) 

By muhiplying the equation (|2.73P j-)-i by (l/'^'')^^i32j (j = 1,2,3), and taking the summation of the 
equations together, we obtain 



\{v'^)-^j:^\B,,A - W"^)-^-vsiEis2,i 



^vii/ "I n/,"i > I r^n^ I ^ 

^ 3 (2-83) 

= iV^-^y ( E Sa.fa, • + E (^21* + AJ21.) • nB,B2i + hit ^ f2,52jB2i) • 

i=i i=i,3 j=i 

It is easy to check that 

\{V'^r''j2B,-hrQ^''\<cjz\B2,? + C\Q^''?- (2.84) 

From the construction of viscous contact wave and (|2.68p . it holds that /^^^ = 0(1)5'^^. Then the terms 
on the right hand side of (12. 83^ can be estimated by 

\{V^'y{l2u + \iT2ix)-riB^B2i\ 

< Ch,T\Bf +C\U^.'\iBf + IB21I') + Ch^rie-'^e-'^ + |(di., ^2., da.)!] (Sf + 1^211') 

< aMB\' + C\U^^m + \B2i\^)- (2.85) 
Similar to ((^^ and (P?55t . we have 

|(y^^)-^(r2it + hhi^) ■ faSaSail < ChMB\^ + CS\U^.^\{BI + \B2i\^), (2.86) 

and 

|(^^^)-^r2i* • E^"2.^2,^2i| < C,,t\B\' + C|[/g| E |i?2,f . (2.87) 
i=i i=i 

Substituting ^^-^^ into ((^^ and choosing N large enough give 

3 \ -, 3 



i(y^-^)-^Eis2,f ) -W^^)-^-^if/SiE 1^2,1^ 



> -C^^tI^P - CIQ^^p - Cd\U^^\Bl - C\U^^\{Bf + S|i). (2.88) 



Multiplying (|2J3)) .^ by F^Ba yields 

Iv^Bl) + h,AV,P)V^Bi - ^NV''-'\V^\Bl 



2 2 ' ' 2 

1 ~ ^ ^ / ^ ^ A ^ (2.89) 



i=l,3 j=l 
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The wave interaction estimations imply 

A3. = X!^^ + X^T + Af^ + (A3. - Af^ - Xg^ - Afl) < Af^ +C\V.\+ Ch^T, (2.90) 

so that 



1 



(2.91) 



The other terms on the right hand side of (12.89^ can be estimated similarly to (|2.8ip - (|2.85p . Then we 
have 

\]-NV^-^VtBl + V^B^h ■ Q^^l < CBl + CIQ^^p, (2.92) 



\V"{ht + \ihx)-fiBiB3\ 

< Ch,T{Bl + Bl) + p\ut \V''Bl + Cp\Ut \Bl + C\%\V''{Bl + Bj), 



(2.93) 



\V'\ht + A3/3.) • ^353'! < Ch,TV'''Bl + C5\U^^\V''Bl + C\V^\V'''Bl 



and 



By choosing /3 and 5o small enough, substituting (|2.9ip - (|2.94p into (|2.89p gives 



2 ^ 



IA3II f.AT j-,2 A3 ,r,>Ar_i|,> I j-,2 



>-C;,,t|S|2-C|Q^^P 

By combining p.82p . p.88p and (|2.96l) and noticing that (5o < 1 and iV > 1, we can get 



W^3)-A'~l|f,S3|^|5^^.|2, 



|A 



>-C„,T|i?P"C,,,T|g^^P. 

Integrating (I^T^T)) over [t, T] x R with i e (/i, T) then yields 



3-^V^Bl-^NV''-'\%\Bl 



[ \B\^dx+ f [ 



\V,\iB-t+Bi) + \U^^m 



dxdt 



<Ch,T / \B\'' +Ch,T e-, yte[h,T]. 



R 



Applying Gronwall inequality to (I2.98P gives 



[ \B\''dx+ f [ \\V^\{Bl+Bl) + \U^^\\B\^]dxdt < Ch,T ei, Vi e [h,T] 
Jr Jt Jr'- ^ 



(2.94) 
(2.95) 



(2.96) 



(2.97) 



(2.98) 



(2.99) 



This completes the proof for the case when k = in Lemma [2.51 The case k = 1,2,3 can be proved 
similarly to the differentiated system, and we omit the details for brevity. □ 
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2.6 Superposition of Waves 

With the above preparation, finally, the approximate superposition wave {V,U,£)(t,x) can be defined by 



where £ = Q 



\u\^ 



Thus, we have 



V 

I it,x) 



V + bi 
Ui + h2i 



{t,x), i = 1,2,3, 



(2.100) 



e = e - ^ u,b2^ + 63 - 



1^2 



(2.101) 



where &2 = (&21, ^22, ^23)* and = bj^- 

From the construction of the contact wave and Lemma |2. II and by noting that cr = e^, we have the 
following relation between the approximate wave pattern (V, U, £ , 6)(<, x) of the Boltzmann equation and 
the inviscid superposition wave pattern {V , U, £, 0)(i, x) to the Euler equations 

|(y,c/,£:,e)(i,x)-(y,c/,f,e)(t,x)| 



< C [I (y«i , C/«i , , ) {t, x) - {v'-' , u'-' , E'-' , ){t,x)\ + \ (rfi , d2 , ds) (i, x) I 
+ (y^^,U^'^,£^'^,@^^){t,x) - (f=^u^^^^^r'^)(t,a;)| 

+ 1 {V'^' , C/^^ ) (t , a;) - {v^^ , u^^ , E^^^ , ) (i , x) I + I (61 , &21 , &22 , &23 M){t,x)\ 



< Ch,T 

< Ch,T 



-(crln(l + <) + (T|lncr|) + - +^C'Dg-7f5^ +,5^3g- 



4^3 



Moreover, the approximate wave pattern (V, [/, 0)(i, x) satisfies 



0, 



f^.t = £( 



V 



£t + (PC/l):r = e( y )X + g£( y + 2^ ^ { 



V 



where P — p{V, Q) and 
Oi = 



Q4 



= — £ 



P - P - (P^bi + Pu ■ b2 + Psbs) 

Qu + Q12, 



4 

— -£ 



A*(e)c/i. A/(e)c/: 



la; 



3 L V 



V 



2,3, 



PC/i - PC7i - ((PC7i).6i + (Pt7i)„ • 62 + {PU^)Eh) 

r,K(e)e, K(e)e,, , 4,^(e)[/i(7i, K^WiUix 



V 

3 



y -) + 3( V 



V 



1=2 
?41 + Q42- 



V 



V 



Straightforward calculation shows that 



(Qii,Q4i) = 0(1)|6|2 



(2.102) 



(2.103) 



(2.104) 



(2.105) 
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3 Proof of Main Result 



With the above preparation, we will give the proof of the main theorem as follows. For this, we will 
first reformulate the problem in the following subsection. The energy estimates will then be given in the 
second subsection. 



3.1 Reformulation of the Problem 

We now reformulate the system by introducing a scaling for the independent variables. Set 

X t 
y = T = -. 

e e 



(3.1) 



In the following, we will also use the notations {v, u, 9){t, y), G(r, y, ni(r, y,£) and {V, U, 9)(t, y), etc., 
in the scaled independent variables. Set the perturbation around the superposition wave {V,U,Q){T,y) 
by 

(</., V, C)(t, y) = {v-V,u-U,E-£,0- e)(T, y), 

G(t, y, = G(t, y, - (r, y, 0, (3.2) 

f{r,y,0 = f{r,y,0-F''{r,y,0- 

Under this scaling, the hydrodynamic limit problem is reduced to a time asymptotic stability problem 

for the Boltzmann equation. 

In particular, we can choose the initial value as 



(<^, V,a;)(T = Ky) = (0,0,0), G(r = = 0. 



(3.3) 



Introduce the anti-derivative variables 



{^,^,W){T,y)= r {^,ij,u>){T,y')dy'. 

J —OO 



Then ($, ^, W){t, y) satisfies that 



^-ir + {P-P) 



V V 



?^(ni-nfi^-nf^)de-Qi-gi, 



\ V 



- I ae.(ni-nfi^-nf^)dc 



Qi 1 ^ — 2? ^) 



K{e)ey K{Q)Qy\ 4 / ll{9)UlUly _ /X ( 8) [/j [/j^ \ 

V ) 3\ V V ) 



(3.4) 



lj,{9)uiUiy ii{Q)UiUi, 



i=2 ^ 



V 



■)- J 6^(ni - nfj^ - nf^)d^ - Q4 - Q4. 



To precisely capture the dissiaption of heat conduction, we introduce another variable related to the 
absolute temperature 

3 

W = W -U ■'^ = W -Y, Ui^h 



then 



1^ |2 

C = Wy-C-^-Uyn 



(3.5) 
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Linearizing the system p.4p around the approximate wave pattern (V, U, 6)(t, y) imphes that 



where 



and 



r $r-*ly=0, 

Z ^ 2 



y2 



3 V 



T-yy 



7^ 



eie.(ni - nfi^ - nf-^)de + a^. - - Q„ i = 2, a, 



V 



'-^.y + Ur-^-^{Uy^)y + ^ey% 



i=l i=l z=l 



3 V ' 



N, = 



4 

-,U.y 



o(i) l$.P + l*.P + ICP + l*i..P 



(C 



N,, 



= 0(1) + + + 



1/2 

z = 2,3, 



Na 



-(p - P)*!,, - ^vl,^ . xp^^ + (^ _ -^Ky 



V 

4 Al{9)Uly _ H{e)Uly 

3^ V V 



V 

fi{0)uiy /i(9)f7, 



ly 



+ ey 

0(1) 



«(0) K{e) , «(e)^ ^'(6),, , \^y\ 



V V 

4=2 

'(-(SlA liT, 12 



-(C- 



. V y ' y2 y ' 2 

l'i>.l' + l*.l' + lCl' + l*..l' + IC.I'^ 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



(3.10) 



We now derive the equation for the non- fluid component G(r, y,^) in the scaled independent variables. 
From (jl.lSp . we have 



where 



Let 



and 



- LmG = -iG, - -Pi(eiG,) - [-Pi(eiM,) - _pf3(aM^^) 
+2g(G,G^^)+Q(G,G) + Ji, 



Ji = (Lm — L]viS3)G'^'' 



JPi(^iG,'-')-^P?(^iGf) 



^""'(^'^'^^ ^ ^Lj^'{PiKi(^^ef +e-c/f )M]}, 



Gi(t, y, = G(t, y, C) - G«^ (r, j/, ^ " G^'^(r, 0, 



(3.11) 

(3.12) 

(3.13) 
(3.14) 
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where G'^^{T,y,^) is defined in IJ^TIU^ 
Then Gi(r, y,^) satisfies 



Gir — LmGi 



Pi(aGj,) + 2Q{G, G^-') + g(G, G) + Ji + J2 - Gf ^ - G 



CD 



(3.15) 



with 



J2 =- 



/ If _ f/Sa |2 



(3.16) 



Notice that in (jXTl) and (|XT5)) . G^i and G'^^ are subtracted from G when carrying out the lower 

T 

order energy estimates because / ||(6y \ Uiy)\\'j^2(^^y)dT is uniformly bounded with respect to e, while 



Both do not give any decay with respect to Knudsen 



|9y ■^||^2(^j^)dT is only of the order of e 



number e in the above integrals. 

From (|1.15p and the scaling transformation p.ip . we have 



fr 



^fy+^-^fy-QifJ)- 



Thus, we have the equation for / defined in 



fr 



Ml 

V 



fy + ^fy = LmG + Q(G, G) + Jf, 



(3.17) 



(3.18) 



with 



= " W^^^y' Vsl^^'^y' + 2Q(M- M^-\G^3) + 2Q(G,G^-^). (3.19) 

The estimation on the fluid and non-fluid components governed by the above equations will be given in 
the next subsection. In the following, we will state the main estimate we want to obtain and also give 
the a priori estimate. 

Note that to prove the main theorem in this paper, it is sufficient to prove the following theorem 
on the Boltzmann equation (I3.17P in the scaled independent variables based on the construction of the 
approximate wave pattern. 

Theorem 3.1 There exist a small positive constants 61 and a global Maxwellian M.^, = M[^,^.„^_g^] such 
that if the wave strength S satisfies S < Si, then on the time interval [-, — ] for any < h < T , there is a 
positive constant ei{6, h, T). If the Knudsen number e < ei, then the problem {S.ll^ admits a family of 
smooth solution /^'''(t, y,^) satisfying 



sup sup 11/^' (T,y,^) - M[y,[/,0](T,?/, 011^2 



re 



< C£5. 



(3.20) 



Consider the reformulated system <\i.Q\ and (|3.15p . Since the local existence of solution to p.6p and 
p.lSp is known, cf. fTS] and [55], to prove the existence on the time interval [-,—], we only need to close 
the following a priori estimate by the continuity argument. Set 



sup ■ 

'-<t'<t 



|($,vI/,W/)(r',.)ir + li(0,V^,C)(r',-)ll? + 



-d^dy 



E 

|a'| = l 



|q|=2 



e e 



(3.21) 



where 9", 9" denote the derivatives with respect to y and r, and M^, is a global Maxwellian to be chosen. 
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Remark 5 In the paper, we simply choose the initial data for the Boltzmann equation i3.n\ l as 



so that 



e ^ ' ' ' e £ 



(3.22) 



(3.23) 



In this case, the functional measuring the perturbation A/'(t) is smaller at the initial time t = j than the 
estimate given in Theorem \3.1\ in the whole time interval when £ is small. 



Note that the a priori assumption (|3.21l) imphes that 



||(a>,vi/,M-)l|2 +||(0,^,C,c.)||2 <c•x^ 



and 



\Giy\' 



d£,dy 



< C 



d^dy 



d^dy 



<CX[X+ I! {Vy, Uy,ey)iQ^- , U , 8^^, t/y^ ^ ) | | L ^ ( ) 



Furthermore, for |a'| — 1 



19" G| 



-d£^dy 



Gy 



-d^dy 



<Cx 



|9"7.P + |9"'(M-M^^^^^^ 



J2 ||9"(^;-^^^^w-C/^^ 61 -9^^)11 



|a|=2 



+ i|9"'(w- v^^^w-c/^^6l-e^^)•a"'(F^^c^^^0^')ll 

|a'| = l 



^ lla"'(?;-y^^w-c/^^ 6* -9^^)11 

|a'| = l 



2 



< Ch,TX^ 



From ([TTT7)) and (|2.103l) . we have 

(t>r - V'ly = 0, 



-Qly - Qly - J £,lGydS,, 

i^^r = -i'-^U,y)y + J a^nfl^dC - Q,y -Q^y- I ^l^^Gyd^, i = 2,3 
0.. + {PU, PU,)y = -[['^Qy)y _ (^^6,^^), 



i'-^U.U.y)y 



V 



~Qiy-Qiv-^ / m'Gyd^. 



(3.24) 



(3.25) 



(3.26) 



(3.27) 
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And from (|2.62p . we get 

Q^y = e[P - - P^^ ~ P^^ - (p^di +p^-^d2 +pfd3)l 

^ 1/ y«i yC-D ySs >^ 

= 0{l)e[\{Vj^\U^\£^\d,^)i\d,\'',V^'~v,,U'^'~u,,£^'-E,)\ + \d,d,,\ 
+ I (C/fl , L/f ^ ^ ) (y «^ - i;. , ) I + I (t/f ^ ^ ) (C/i^^ , , di, , ) I 

■— Ql3 + Ql4, 

where Q13 represents the wave interaction satisfying 

Qi3 ^ 0{l)e[\{Vj'\U^\£^\d,,){V''' ~ v,,U''' ~ u,,£''' ~ E,)\ 
+ 0{l)e'[\iU^,\U.^-V,^-){d,,V^^~v^)\ + \{U^^ 

+ |(f/f^K''^)(t/i^^K''^rfl.,d2.)|] < c„,t e-^e-^, 

and Qi4 represents the terms related to the hyperbohc waves di (i = 1, 2, 3) satisfying 

Qi4 = 0(l)e [1 (y/i , [/f 1 , , d,.) I |2 + I 

+ 0(l)£2[|(d2.,,di,d2.,V/i(i2x,t/f^(ii,)| . 

Then we have 

' \Qiy?dydT <C I [{\Qi3\^ + \Qi4\'')dydT<Ch,Tei. 



Similar estimates hold for Qiy [i = 2,3,4). 
By (|2.104l) . straightforward calcuation gives 

Qiy = e[P - P - [p^hi +pu ■ &2 +Pe03)\^ - -e {- 



V 



V 



0(l)e^ [|6.n(t4,t/i.,4)| + |6.||6. 

1=1 
3 



Hence, 



Similar estimates hold for Qiy {i — 2, 3, 4). 
Thus from the system (|3.27p , we have 



|l(</),,Vr,C^r)f<C„,Tx', 



and 



(3.28) 



(3.29) 



(3.30) 



(3.31) 



(3.32) 



(3.33) 
(3.34) 

|i(0,,Vr,Cr)||' <C||(0,,Vr,C^r,t/r-V^)|P <G>,TX'- (3-35) 

Now we want to obtain the estimates on ||9"(0, V', C)IP for |a| — 2. For brevity, we only calculate 
\dyy{v — V^^,u — — 0"^^)|, and the others can be estimated similarly. From (|1.3p and (I2.5ip . we 

have 

I dyy [v -V'^'',u-U^\6 

= o(mVy%^u^l,Q'yl, (y/3)^ {u^-f, (6^^)2)1 • \{v ~~ v^\u~ u'-\e~ e^^)| 

+ 0(1) I (Kf ^ , U^-^ ^Ql^)\.\{Vy- If ^ , Uy - U^^ , Oy - 9 f ' ) | 



+ 0{l)\{Vy~V.^\Uy 



U^'.ey-Ql')? + 0(1)^2 J \v^iO\-\dyym- (3.36) 
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Therefore, we have 

+ ||a"(F^^,f/^^,e^^)|p + ||a"(di,d2,d3)f 



621, &22,&23, ^3)11 



c 



(3.37) 



Finahy, by noticing the fact that / = M + G and i^-^^ ^ ^. qSs^ |337| . with |q;| = 2 yields, 



a<~i |2 



-d£.dy < C 



|9"(M- M^3) 



-d^dy < Ch,Tx\ 



(3.38) 



where in the last inequality we have used a similar argument used for p.26p . 

Before closing the a priori estimate (I3.2ip . we list some basic lemmas based on the celebrated H- 
theorem for later use. The first lemma is from 1141. 



Lemma 3.2 There exists a positive constant C such that 



M 



-d^ < C 



M 



M 



P 



M 



M 



where M can he any Maxwellian so that the above integrals are well-defined. 

Based on Lemma \3l2\ the following three lemmas are taken from [29]. And their proofs are straight- 
forward by using Cauchy inequality. 

Lemma 3.3 // 9/2 < 6^, < 0, then there exist two positive constants cf — cf(v,u,9; 
Vi,,Ui,,0-):) and 770 = ?7o(w, w, 0; w*, w^,, 0^) such that if \v — + \u ~ u*| + \0 — 6'*| < rjo, we have for 



d£,>o- 



Lemma 3.4 Under the assumptions in Lemma \3.3[ we have for each g{£_) G 



Kiel) 
M 



M 



d^, and 



Kiel) 



de 



Lemma 3.5 Under the conditions in Lemma \3.3\ for any positive constants k and \, it holds that 



giPi(iei''g2) 



dc 



gl|CI'g2 



dCI < Ck,x 



A|giP + A-i|g2p 



d^. 



where the constant Ck.x depends on k and A. 



3.2 Energy Estimates 

To close the a priori estimate (I3.2ip and to prove Theorem 13. 1[ we need the following energy estimates 
given in Propositions 13.11 and Proposition 13.21 First, the lower order estimates to the system p.6p and 
p.lSp are given in the following Proposition. 

Proposition 3.1 Under the assumptions of Theorem \3.1l there exist positive constants C and Ch,T 
independent of e such that 

|GiP 



sup 

!}<T1<T 



■{Ti,y,C)didy 



\Uf^\{:i',W)f + \\{%,^y,WyX,^r,WrW 

<Ch.,Te f \\{^,W)\\'dT + C I ||a"'(0,^,C)|pdT 

•^■F \a'\ = lE 



dT 



Gi pd^dydT 



^ II 



|a'| = 
^(l^l)iaa'P^|2 



a" G\''d£,dydT + Ch,T e^. 
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For brevity of the presentation, we also put the proof of Proposition 13. II to the Appendices. 
Then we perform the higher order estimates. Firstly, we apply dy to the system (I3.6P to get the 
following system for (0, C)- 



i'ly = 0, 



-A^5 - Qly - Qly, 



3 



3 „ 3 

+ a6(ni - nfi^ - nf%d^ + Ns- Qiy -Qiy + Y. u^{Q,y + Q,y), 



(3.39) 



where Z is defined in p.7p . Here, the linear terms are 

Hi : 



- {y)y<f> + (^).C - (^^t^l«),C, 



(3.40) 



i?, = — — — r^C^^aCa - ( — — ly't' - \—{^^''v)y^^ * = 2, 3, 



y2 



V 



(3.41) 



if4 



3 



V(e)c^^^4Me) 



.^(6) 



V 



V 



(3.42) 



and the nonlinear terms are 



N. = 



P~p, , 4 u(0) /i(e) 



/M^?) M(e)_^M0), m'(©) 



+ 3^— -^)^^^+3^^H— + 



y " 3' 1; 
= o(i) [|(0, V, OP + l(0y>i„ Cy)P + IV'iyyP 



1/2 



(3.43) 



i+4 — 



M(0) 



o(i)f|((/),^,C)P + !(</-„ V-.y, OP + 



, 1 = 2,3, 



\ V 



(3.44) 



and 



-{p- P)ijiy + -^{p- P)(j) + 2[— 



y2 ^ 

k(6I) k{Q) 



i=2 



){-^Uly^ly + ^U^y'lp^y^ 
4=2 



0(i)[|(<^,V',C)P + l('^y,V'y,C,)P + ICy,P 



V 



y y2 



(3.45) 
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To derive the estimate on the higher order derivatives, applying dy to the system (|3.39p . gives 



(pyr - tPlyy = 0, 

Z 2 
Vlyr - y4>yy + ^Cra + 



.4 Me) 



3 V 

CD ttSs 



C,yr + Zipiyy + 

3 



^l(ni - Hii - Ili^)yyd^ + A^Sy - Qlyy ~ Qlyy, 

"("(-^2+4)?/ Qiyy Qiyyi ^ — 2, 3, 



yy 



^,\^{u,-n?,^-uf%yd^ 



"'— '' ^ t— 1 

3 

(94^^ — Qiyy + E I^C^i 



-TVs 



(3.46) 



where 



and 



l^U.yCyy ~ K^)y^y - {^)yy^ + ( ^ ) , 

2 .8M'(e) .4^(9) 



i+4 



/i(e);7,y, fi'ie) 



V 



^ty'^yy ^ ^\ y2 ) yVy 



y2 I'VV 

-^\—;^^^v)yCy-[—^U^y)yyC. « = 2,3 



y2 /yy^ 



8 — Zy^ly 

3 



V 



2. V(9)[/i, 4^.(6) , .1 



i=2 



V 



3 

y V 



(3.47) 



(3.48) 



(3.49) 



By using the above two systems and the equation for the non-fluid component, we can reach the 
following proposition for the higher order energy estimates. 



Proposition 3.2 Under the assumptions of Theorem \S.ll there exist positive constants C and Ch,T 
independent of e such that 



sup 

T<-ri<r 



b,^P,C^y,^yXy){ri,W+ E / [ K^i^i^V^^d^dy + ( j \^lL{r^,y^^)didy 



Kiel) 



~ l<|a|<2 

< C{5 + Ch,TX) 



l<|a|<2' 



'^\d-G\Hidydr 



\GiYdidydT + C{5 + Ch,Tx) / \\{<l),^kX)rdT + Ch,T e 



Again, the proof of Proposition 13.21 will be given in the Appendices. 

By combining the above lower and higher order estimates given in Propositions 13.11 and 13.21 and 
choosing the wave strength 5, the bound on the a priori estimate x E^nd the Knudsen number e to be 
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suitably small, we obtain 

-^W+r[ E \\d-ic^,i^,Or + \\^/pS'\{^,W)r]dr+ r f f d^dydr 

•'i 0<|q|<2 ^ Jh J J * 

l<|a|<2 ~ * 

Therefore, we close the a priori assumption (I3.2ip and then complete the proof of Theorem 13. II 

4 Appendices 

As mentioned before, since the proofs of Propositions 3.1 and 3.2 are technical and long, we put them in 
the following three subsections. 

4.1 Proof of Proposition 13.11 

Proof of Proposition 3.1. Firstly, from the fact that 

< V^-^ < V*, V'^^ e (min{'y*,w*},max{t;*,w*}) and v* < V'^^ < v+, 

we have 
Thus 

^ <V < 2v+, if e « 1 and S'^" < 1. (4.1.1) 
From (|2.57P ;^ and ()2.57P o. we can obtain 



where ai = p+ + s§f + = p* + s'^v* 
Since 



3 



ia^-slV''^)- Jei^f'd^, (4.1.2) 



< p+ < (ai - slV'^") <p*, and \J ^.fnf'd^l < C(5^^ (4.1.3) 
we have 

p+ - CS^' < Z^' <p* + CS^', 

and then 

Y <Z^^ < 2p*, if (5^^' < 1. (4.1.4) 
Now we estimate Z defined in p.7p as follows. 

4Me«0<^,,ci. 4M(e^^)£/f,^ 



z = z^-^ + (p^i -p,) - + (p' 



_ pRi pCD pS3 



3 T/fli 'V ^ ^ 3 yCD 



4 



M(e)c/i, A^(e«o< M(e^^)c/f,^ ^(e^^c/f^^ 



y yCD ySi 



where is related to the hyperbolic waves and the wave interaction terms given by 

3 

Q, =0{l)Y,\{d,,h,)\+0{md2yMy)\+0{l)e~^ 

= 01 -+£4+ +e - , 
cr (7^ 



(4.1.5) 



(4.1.6) 
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with (T = e 5 . 

From the properties of the approximate rarefaction wave, we have 

0<P«i _p, <p__p,, and -^^^—2-2lL < Ch,T 



And the properties of the viscous contact wave imply that 
pCD _p* < Ch,T e, and ^ 



3 yCD - ^'^'^ 



< Ch.T 



(4.1.7) 



(4.1.8) 



Substituting (14.1.21) . (|4.1.6p . (|4.1.7p and (14.1.81) into (|4.1.5p . we know that there exist positive constants 
c and C such that 

0<c<Z<C, (4.1.9) 
provided e ^ 1 and the wave strength (5 <C 1. 

Step 1. Estimation on || ($, VK)(t, OIP- 



By multiplying dHli by ([321)2 by ([SSIa by y'^^^'i, by gfyW^ respectively and adding 

all of them together, we have 

/l($, W^), + /2(*, W^) + /3(«'y, Wy) ^ /4(*, W^, M^y) 

-(A^i - Oi - Oi) + 1^^^ ~ " Q')*^ 

'i=2 

3 3 



3Z2 



(4.1.10) 



where 



/!($, ^, M^) = - + + i_ ^ ^2 ^ 



j=2 



^3Z2' 



3Z 



2k{Q) 

3Z2y 



(4.1.11) 



2 2M^ 

/4(4', a>„ 4-,, M/,) = - ^w-vi;i([/i, + z,) - — E c/..*. 



3Z2 



2*1 u(e)y^^ u'(e)T/^^ , j A 



+ 7: 



3Z 

4 M(e) 

3^ Z 
4a^ 
3 Z 



y2 



3 



i=2 



2k{Q) 



W^Z^ 



WZ'^ 
11 



W{Uy ■ *), 



3V^Z2 

and Ki denotes the non-fluid parts given by 



4j 



1=1 



Ki = / ei(ni - nfi^ - nf-')de - / ^i^^(ni - nfi^ - nf^)de 



214^ 
"3Z2 



(4.1.12) 



(4.1.13) 



37 



From (|4.1.9|) . we have 



c($, Wf <h< C($, Wf, 



for some positive constants c and C . 

Now we estimate /i, [i — 2, 3, 4) term by term. Note that 



S3 j C?nf^dC, 



and 



Zr 



where Z^^ = P^^ - - 
3 



= P 



CD 



3 yc^^ 



and 



i^n{e)Uiy M(e^^)f/it M(e^^)f/f,^ M(e^^)f/f,^ 

3 



yCD 



2/6^ ef^ e|^\ 2/8]^ e^^T//^ qcj^cj e^^F/^^ 



3\v yc'^ ^■33/ 3V y2 (v^^)^ (t/cd)2 (^53^2 

4f^^{Q)u,yK M(e«0<y/^ M(e'^^)c/g^K'^^ M(e^^)t/i;^F/- 



T/2 (yfii)2 

'3 1 y 



yCD 



VS3 



Then we get for the first term of I2 that 

2 1 1 V 

3Z^'y~2^Z^^ = 6Z^'y+2Z^^^ 



+ ^ (C^f; + t^S" + ^2. + 621,) + ^ (z.^- + + E 

Z— 1 



6Z 

6Z2 



where 



Qr = ^[{Z- _ 3,2(^ _ ^S3)] + + + rf^^ + j,^^^^ 



6Z2 



Firstly, straightforward calculation gives 

-(ai - sjV'^') + 3s2F^3 > _ ^jS3 > 3p^^ jf j53 ^ 



I J ei^f'd^l < and I J Cinf^dCI < CJ^^|(7 



S3 1 

ly I 



And (|4.1.18l) implies that 



2^Z' 
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provided i5 <C 1. 

By the definition of the approximate wave pattern defined in (|2.100p and (|2.10ip . we have 

Qei = 0(l)e[|(ef \ - i;,, (7f ^ - til,, T/^^' - 1;% C/f^ - uDI + |(ef Ml, 

+ I (ef ^, C/g^) 1 1 - V, , f/f 1 -u,,,V'^-v*, C/f ^ -ul,d,,d2)\ 
+ I (ef (7it ) 1 1 - , t/f ^ - i/i* , - t;* , (7f ^ - ul , di , d2 ) I 

+ \id2ud3t)\ + |5i||(ef\ ef^,ef^)|] + ^(^3 -tj.b,- M). 

< Ch,T e + C|(fe2,62r,63r)|• 
On the other hand, direct calculation gives 



i63 



+\ugf 1 1 (1^^^ - V, , c/f 1 - t/i, , - 0, , - V* , f/f ^ - , e^^ - r , di , d2, rfa) I 
+ 1 c/g* 1 1 (1^^^ - 1^* , L/f ^ - lii* , e«^^ - 0„ T/^^ - 1;* , [/f ^ - lit , e^^ - r , di , d2 , d3 ) 



+ |C/ll^,||(di,d2,d3)| + |(6l,62,&3)||(C/S,C/!.r,t^S)l + \d2.t 
< Ch,T e + C\b21yr\- 

Similar estimates hold for Q62, Qei and Qez- 
Moreover, we have 



3V 



321J/T 



and 



3 

Now from (|iX^ . (|4.1.19p and (|4.1.2ip - (l4.1.24p . we have 



-)r\ < Ch,T e\ 



J \Q7\^\dydT < Ckt e j J 'fldydr + C J \{b,by,br,byr)\\'i/i\^dydT 
< Ch,T e j^^ J ^IdydT + C \\{b, 
<Ch,T£ I / 'ifldydr + Ch.T I H^-i || ^ H^-iyH 



< /3 \\'^iyfdT + Ch,T^p e 
Substituting (|4. 1.251) into (|4.1.20p gives that 



"^{dydT. 



I \U^y'\^ldydT - (3 \\^,yfdT-Ch,T,pej^^ j ^Idydr, 



provided 6^^ <^ 1. 

On the other hand, we have 



Similar to (j4.1.26p . we have 



*2 



1=2 



i=2 



dydT = I I —ZrW^dydT 



/V si ( ~ ^'C^il; + S3 1 eiTlf^di + + Z^- + Qe) W^dydT 



\U^y'\W^dydT~l3 



\WyfdT-Ci 



W^dydT. 
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Thus we complete the estimation for I2 as 



j^^ j h{^, W)dydT > j^^ j |C/§ W)\''dydT 



/5 / miy,Wy)fdT-Ch,T,p e \ \ \{^,W)\^dydT. 



(4.1.27) 



Direct computation yields that 

Now we estimate /4 in (|4.1.12p . Note that 



Uir + Zy = - I ei^?iyd^ - / eXd^ + Qly + Qly 



Thus 



^infi^dC + / C?nf^de - Qly - Qiy 



3^2 



)y{Ql+Ql) + {---)y 



< 



(4.1.28) 



where from now on, /3 is a small positive constant to be determined and is some positive constant 
depending on /3 but independent of h,T and e, while Ch,T,i3 depends on on h,T,l3 but independent of e. 
Similarly, we have 



^4 



E 

1=2 
3 



3Z2 



CiC^nfi^de + / CiC.n^M^ -Q^y-Q 



•■lly<- 



1=2 i=2 



ly^ 



3Z2 



< !(*,„ Wy)\' + Cp5'-^ \u!^m^.W)\^ + Cu,T,p e ^ |(*., W^)|^ + (• • • ), 



1=2 

Next, we can show that 
3 

fy, i^'y^'^y 



i=2 



i=2 



$,,-^^^^4 {y/y+lJy^)\Y,lJ,y^, + {--\ 



(4.1.29) 



where we have used the fact that Ui — Ui + fe^+i = Uf^ + h+i — 0(1)£2^ (i = 2,3) and that Uiy 
Ufy'' + ih+i)y = Oil)si,t = 2,3. 
By Holder inequality, we have 



11 



J2ll<m%,'^y:Wy)\' + C0S'''\U^%'m:W)\^ + Ch,T,p e\{^,W)\' 



i=4 



Hence, we have 



h{-9,W,%,-9y,Wy) 

< /3|($„ Wy)\^ + (Q<5^^ + C,^T e)|f/f,^||(vE', W^)|2 + C;,,t,^ £|(vI/, Ty)|2 + (. . . )^. 



(4.1.30) 



(4.1.31) 



Now we estimate the nonlinear terms ^A^i^i by 



V 



N,^,dydT <CJ^J [\%\^ + |vl/,|2 + |^|2 + i^^^^i^ 
<Cx r||(<i>„*^,C,*i,y)lPdr. 



dydr 



(4.1.32) 
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Similarly, we have 



J, J V''^N,^., + —NidydT<Cx my,^yX,'fyyXvWdr. 



T " 1=2,3 



We now turn to estimate the terms -^^'iQi, - Ya=2 V^^'Qi'^i and -^{Qa - Yfi=i UiQi)- From 



(p:^ and (E^, we have 



V - „ - 2W - - 

i=2 i=l 



4=2 

< CK*, W^) 1 1 (Qn , O2 , Os , Q41 ) I + q (*, W^) 1 1 (Qi2 , ^42) I . 
And from (|2.65p . we have 

\{^,W)\\{Qn,Q2,Q3,Q4i)\dydT<C f /|(*,W^)|e 

R Jr 



(4.1.33) 



' dydr 



<e \\{^,W)fdT + Ch^T e-—. 

On the other hand, from p.62p - p.64p . Lemma 2.3 and noting that tr = e^^ we have 
\{^,W)\\{Qi2,Qi2)\dydT 



(4.1.34) 



R 



dydr 



(4.1.35) 



\L^{dy) + £\\di\\L^(dy) 



dT 



, £ , 3 



<Ch,Ti-)' \\i'f,W)h2idy)dT<e \\{^,W)\\'dT + Ch^T s-^. 
Now we estimate the terms -EL2^'^'Q«^i ^nd -f|i(Q4 - ELi ^^^QO- For this, from 



(PJUi| and dUnSl), we have 
V 



R 



2H/ 

^*igii-^(Q4i-c/iQii) 



dydr 



<C r f \i^i,W)\\{Qii,Qii)\dydT<C Tf \{M'^,W)\J2h\^dydT 

~ ^ e R. i — 1 

<C ||(*i,W^)|U2(rf,)^||fo,||i»||6,|U2(rf^)dr 
<Ch.Tei[ \\{^i,W)\\dT <e [ W)\\^dr + Ch,T e. 



(4.1.36) 



On the other hand, from Lemma 12.51 with the properties of the hyperbolic wave II, we can show that 

?12 - ^^'Q'*' - 15(^42 - U1Q12 - Y '^'Q')] '^^^^ 



R 



z 



R 



3Z 9Z2 



V 



V 



V 



V 



(4.1.37) 



214^/^(9)9^ m(0)0i;\ 8W f ^{e)UiUiy fi{e)UiUiy 



zw / 



3Z2 V V 



V 



aw n 
9Z^V 



V 



V 



2w ^ ^^^{e)u,u,y ^iie)u^u, 

i=2 



21^^ // 
3^2 ^ V 



ly 



V 



V 



dydr Y^hi- 
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For brevity, we only compute /51 in the following because I^i {i = 2, 3, 4, 5) can be estimated similarly. 



/51 



, — Qz^)\-^ —r'^'' 

u/^(0), ^/M(e) M(e)v7 



( 



dydr 



' r f /4V^*i 8UiW\/fj.{e) M(0)\r. 1, , ri r2 



(4.1.38) 



We can further have 



4 



Oil] 



R 



W^,) I + I C/i„ e„ Z,) 1 1 (M/i, W^) Ij dydr 



</3y^ [Il(*i„w^,)lli2(,,) + ||^|c/f,1(vi/i,w^)lli.(,,)_ 

dT 



(4.1.39) 



</3 



ll(*i.,w^.)lli^(<i,) + llVl^i1?K*i'^)lli^ 



Similarly, we have 



iL = 0(1) 



R 



\{'fi,W)\\Uiy\\{h,b2,b3)\dyd7 



<P \\^W!^\{'^i,W)fdT + Ch,Te II (*i, W^)fdT + C,,,T,/3 



(4.1.40) 



By integrating (|4. 1.101) with respect to y and r, then combining all the above estimates, and choosing 
/?, (5'^^, £ and x small enough, we have 



($,*,M/)(r, ■ 



c/f;i(*,w^)ir + ll(vi/„w-,)ll 



dT 



(4.1.41) 

+Ch,TeJ^^ \\i^,W)\\^dr + cJ^^ J Kidydr + Ch^T,0 el 

Now we turn to estimate the microscopic term J Kidydr in (I4.1.4ip . Here, we will only estimate 
Kii - J ^-^i y ^i(ni - nfi^ - Ilf^)d(,dydT because the other terms in J J Kidydr can be 

e e 

estimated similarly. Let M^, be a global Maxwellian with the state (w^,, u^, 6'^,) satisfying ^6 < < 6 
and |w — Wi,| + |u — u*| + 16* — 0^1 < r/o so that Lemma [3.31 holds . and Lemmas 12.31 and 12.41 hold with Mg 
being replaced by M^,. Note that the above choice of the global Maxwellian can be obtained if the 
total wave strength S is small enough. By the definition of ni,nf]^ and nf^ given in (|1.20p . (|2.2ip and 
(j2.58p respectively, we have 



Hi - nfi^ - nf' = L^[Gr - ^Giy + ^Pii^iGiy) - Q(Gi,Gi)] 

-Li^'[2Q(Gi, G^^ + G^^ + G^^)] + J3 + J4 + Js- 



(4.1.42) 
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Here 

= ^^Tc^L- 4Pf-[a(^^^e^- +e • (4.1.43) 



3 .-1 r.,ci.r. .l^-f^'^^P 



and 

J5 = L-i[- ^G«^ + iLj^iPi(6Gf ) ~ Q(G«SG«0]. (4.1.46) 

By (|4.1.42p . we have 

^ii = -/jy'|«'i / e?Lj^^(G.)dCd2/dr + ^^^^ 1 Ci2Lj^i(Gi,)dCd2/dr 

+ jl j efLM'[2Q(Gi, G^- + G^^ + G'^^mdydT - Yll J f.Jsdidydr ^^'^'^'^^ 

~L I^I ^i^^rf^c^yrfr ^ m^llj e^hd^dydT ^ ifj 
For the integral we have 

= - /Y 1*1 / ei^MiGir)didydT 



1=1 



' / / 1*^ / ^'I^M (Gf ^ + G?^)d^dydr =: if + if J?. 



(4.1.4 



Note that the Imearized operator Lj^^ satisfies, for any h £ *Tt^, 

{L^g\ = L^{g,) ~ 2L^{QiL:^g,M,)}, for , ^ T,y. (4.1.49) 

Then we have 

Kll = - / ei(J^^G,UCdydT^2 / e?LMHQ(LMGi, M.)KrfydT 

= - Jj [^^iei'L^{Gi)]iT,y,Od(dy + j'^ j j ei^^iG ^)d^dydT 



-2 / ei^^lQi^^mGuMrmdydT. 



(4.1.50) 



The Holder inequality and Lemma 13.41 vield 



eM(Gi)dep<c/ i^de (4.1.51) 



Moreover, from Lemmas 13. 2113. 4[ we have 

J e2L-i{Q(Lj^iGi,M0K<c(^y' ^|L-MO(L-iGi,M.)}prfe) 



<^( /4^|LM^G,Pdc)'-(/«|M.Prfe)' (4.1-52) 

-1/ 
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Combining (I4.1.50|) - (l4.1.52p gives 



|GiP 



\\i^^r,J\U^:\^i)rdr 



iT,y,^)d^dy + Cp 



Kiel) 



(4.1.53) 



+ Ch,Te \\^l\\'dr + ChAe~' +X) ||(0r>r,Cr)|rdT, 



On the other hand, by p.l3p and (|2.20p . we have 

Kl? = -fJ 1*1 / eil^M^G?' + G?^)d^dydr 

< c m |*i|[|(e«s[/«v, 9^/^,(7,^,^)1 + \ie^\u^;,&^^,u;f'')\\ivr,ur,er)\]dydT (4.1.54) 

<Ch,Te r \\^i\\^dr + C{Ch,T e + S) [ \\i<Pr, ^rXrWdr + Ch,T 



which, together with (j4.1.53p . imply 

</3[||*(r,-)l|N 



-Cb 



||(vl/i„^|C/^;;'|*i)||'dr 
iT,y,C)d^dy + Cp 



Ch,T e = 

Kiel) 



Gi\^ d^dydT 



(4.1.55) 



+Cn,Te ||*i|prfr + C[5 + C„,T(£' +X)] / ||(0r,V'r,Cr)||'dr. 



We now turn to Kf^. By (I4.1.49p . we have 



K'' — - 
-'^11 — 



< C 



m(^^^^)yj e{^^miGi)d^dydT + 2 j^j ^^ij^ j ei^m{Q{^mGiMv)}didydT 



\^^y\ + \^^\\{Vy,Zy)\ + \^^\\{Vy,Uy,ey)\]{^ 

Kiel) 



'^-Hiei) 



Gi\^dA'dydT 



</3 \\{^^y,j\u,%^\^^)rdT + c^ 



Gi\^ d^dydT 



+ Ch,Te \\^>irdT + ChAe-' +X) Wi^y.^yXyWdT 



To estimate Kf^, notice that 



Pi(eiGi,) = [Pi(aGi)], + ^ < eiGi,x, > Pi(x,,)- 

Then from (|4. 1.491) and (|4.1.57p . we have 



(4.1.56) 
(4.1.57) 



^11 — 



{^)v I ei^^[Pi{^iG,mdydT 



"2 J e?LMHQ(LM[Pi(aGi)],M,)}dCd2/dr 

</3 j{\U^^\\^^\^ + \^^y\^)dydT + Cp J^^ J 

+Ch,TS J \-f,\^dydT + Ch.T{e'^ +x) J 



(4.1.58) 



Kiel) 



Gi\^ d^dydT 



Mv.'>PvXy)?dydT. 
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On the other hand, 



Kti < C 



< C 



!^\L^{QiG,,G^)}\'d^ydydT 



1*1 



Gi\^ d^dydT < Cx 



Kiel) 



(4.1.59) 



Gi\^ d^dydr, 



and 



i^fi < C 



1*1 



Kiel) 



</3 



\U'^\\^i\'dydr + Cp 



Kiel) 
Kiel) 



G^-'pde)' + |(ef\e^^,t/«;, {7^^-13)11 ^^^^ 



Gi pdfdydr + Ch,T s 



l^ipdydr. 



By (|4.1.44p . we have 



(4.1.60) 



K° — - 



^ U^^^ 1 1 

- ( — ^M - ■jT^L^^Sa)^^^ + (-Lm Pi - -^Lj^sPf')(eiG^') 

5 

- L^[2Q{G^^,G"' + G'^^)]}dedydT := ^ iff 



(4.1.61) 



'^11- 



For brevity, we will only consider the term Kfl because similar estimates hold for the terms Kfl (i 
2, 3, 4, 5). Direct calculation yields 



K 
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< c 



z 

l*ll( 



e?Lj^i [2Q (M^^ - M, L-is3 (Gf ^ ))]d^ 



t^(|g|)|M^-^ -Mp 
Ml 



Kiei)|Lr,is3(Gf^)p 



< 



C/ / |*l||(^'-F^^^i-^7^^6'-e^^)|( 



^-^(leDiGf^p 



d^Ydydi 



d^)^dyd7 



<CSJ^ J \U^^\\^i\^dydT + Cd j \{%,-^yX)?dydT + Ch,T e J \^i\^dydT + Ch^t . 

(4.1.62) 



By (|4.1.45p . we have 



.G 



CD\ 



'V yCD M.'^'^^y J 

[L^fQ(G^^, G^^) - L^ie.g(G^^, G^^)] 
+ [iLj^iPi(6G^^) - -LjL^^c.Pf^(eiG;;^)] - L^mG^\G^^)]]didydr := ^if 



For illustration, we only consider KH as follows. 



< C 

< Ch,T e 



Z 



e? 



CD 



)Li\;r G 



TjCD 

7 -L'tv 



/C'UU _ (~tCD 



yCD J M. '-^y 

dS^dydr 



jCD 



pi CD 

^CDi^y 



l*i||(w-i^^^,u-c/^^,0-e 



|(e;;,^,[/,^,^)| + |(e^^,[/,^^)||(^;„u„0,)| 



Ivfil^dydr + Ce / \\{^y,^!y,OrdT + Ch^tx \\{4>v,i^vXy)rdT + Ch,Te 



dydr 



(4.1.63) 
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Then, by (|4.1.46l) . we have 



< C 



dydr 



(4.1.64) 



< Ch,T e J l^ipdydr + Ch^t e 
By collecting all the above estimates, we have 

||v|/(r,.)f + r \\{^,r.\^\^l)fdT 



+Ch.T 



Kiel) 



Gi\^d^dydT + C{x + P + S) / my,^yXWdT 



(4.1.65) 



: / W-^ifdT + CiS + Ch^Tx] J2 L \\d"\(b,^jX)fdT + Ch,T 



Therefore, we have 

||($,vI,,M^)(r,.)|P 
|GiP 



|c/f,1(*,w^)lP + II(*.,W^.)I 



dT 



< C 



iT,y,Od^dy + Ch,Te IK*, W^)||^dr + C/3 / \\{^r,Wr)rdT 



+Cix + P + S) / ||($„C)irdT + C/3 



r 

+Ci6 + CH,TX) E / 



Kiel) 



Gif d^dydT 



(4.1.66) 



lia" (0,v,C)ll'dT + a,T,^ 



where we have used the smallness of S, /3, e and x- 
Step 2. Estimation on ||$y(T, •)|p. 

Note that the dissipation term docs not contain the term |l$j^p. To complete the lower order energy 
estimate, we have to estimate From 



4 M(e) 

3 V 



V 



3V 



2, we have 
2 



3 V 

Muhiplying (j4.1.67p by $y yields 



Uly{Wy + C/y • *) - A^l + gi + Ql + / ^li^l - H^^^ - 0^^)^^ 



(4.1.67) 



.2^(0). 



Z 



4/x'(e) 



3 y 

where we have used the fact 



Uiy{Wy + Uy-^)-Ni + Qi + Qi+ / (Hi - n^i^ - n^^)de $„ + (•••) 



r 2 
3F~ 



31/ 



-.W„ 



(4.1.68) 



Integrating (|4.1.68p with respect to y and r and using Holder inequality give 



2/i(e) 



3F 
< Ch.T e 



(T,y)dt/+ / / —^Idydr 



\m?dT + C / \\{-^y,WyWdT + C 



(4.1.69) 



-C5 
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By (Umi, (1X51) and (I2.104p and the Cauchy inequality, one has 

{N! + Ql + QDdydr < Cx IK*,, *„C)lPrfT + Ch,TX Uiyfdr + Ch,T (4.1.70) 
Now we estimate the last term on the right hand side of (|4.1.69p . By (|4.1.42p . we have 



< c 
+ c 

+ c 



e,J^^{Gr)dC\'dydT + C I III ei^Mi^Giv)d^\'dydT 



e,L^[-P,{^,G,y)]dedydT + C 



eiL^[Q{Gi,G,)]dedydT 



+ 



C jl j J C?(|J3p + |J4p + |J5pKrfydT:=^K^. 



Then we can obtain 



<C 



< C 



Similarly, 



< C 



< C 

< C 



^-KmGiy? 



M 



dS^dydr 



-d^dydr. 



'^"(l^l)l^^l -d^dydr + c I ||(e«;,<„e^,^,t/,^,^)|pdr 



'-d^dydr + C fj I ^l^dGf p + IG^^pRdydr 



< C 



^-\\mG . 



^' d^dydr + Ce I \\{(j)y,iljy,Cy)fdT + Ch,Tei 



Moreover, 



< C 



< C 

< C 



'ly) 



M 



l2v^,2u^,.20^,] 



/ '^'\\mtM^2v.^2u^^2e^]d^\dydT 



-i(|e|)|iPi(eiGi,)P 



HmG '2 



d^dydr < C 



<mGiy\' 



[2vt,2ut,2et] 
12 

^' d^dydT + Ce I \\{(t>y,ipy,Cy)fdT + Ch,Tei. 



■d^dydr 



From Lemma l3.2| we have 

^-Hie|)|Q(Gi,Gi)p 



< Ch,TX 



d£,dydT < C 



(4.1.71) 



(4.1.72) 



(4.1.73) 



(4.1.74) 



■dS^dydr 



-dS,dydT. 



(4.1.75) 
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Using similar idea as for (I4.1.60p and (|4.1.62p . we can obtain 



and 



<Ci6 + e) / ||($y,*y,C)||'dT + C,,,T ei 



-C((5 

I h I I ivi , 

|a'| = 

Thus combining (I4.1.69p . (14.1. 70p and (|4.1.78p yields 



+C{S + x) j[\\{<Py,i^yXy)\?dT + C mi '^\d"G\'d^dyd7 



+CS / ' ||y^*fdr + C[<5 + C,,TX] /' / / ^\G,\'d^dydT. 



Step 3. Estimation on the non-fluid component. 



(4.1.77) 



Substituting the estimations of K2 {i — 1,2, ■ ■ ■ ,6) into (|4.1.7ip . we obtain 

<C{5 + e)j^ \\{<i>y,^y,C,4>y,i'y,Cy)\?dT + Ch^T 

e |Q!'I = 1 e 



(4.1.79) 



The microscopic component Gi can be estimated by using the equation (13. 15^ . Multiplying p.lSp by 
tIt gives 



+g(G,G) + Ji + J2-Gfi -G?^|^ (4.1.80) 
I - Pi(^iG„) + 2vQ{G, G^') + vQ{G, G) + vJi + VJ2 
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The Cauchy inequality implies 



< 



32 

+C 



- |q|=2 

{Vr,Ur,9r)\^ + \ {Vy , Uy , ey)\^ 



dydr 



< 



32 / h. 



m J '^^^M^' d^dydT + CiS + C^^Tx) jl E W'MXWdr + Ch,Te-^. 



|a'| = l 



Notice that Pi(6Gy) = - ELo(Ci Gy, Xj)Xj- Then we have 



e i=0 E 



< 



Kiei)Gf 



d^dydr + C 



32 Jh J J Jh J J 

+C(5 + C,.,TX) / \\dy{<p,4,,CWdT + Ch^Te"^, 



d^dydr 



where we have used the fact that 



|(6G„x,)P< 



-dt 



Lemma l3.2l -Lemnia 13.51 and wave interaction estimates imply 



Q{G,G)^d(dydT 



< C 



d^ 



< 



< 



32 

a 
32 



/ / '^^^d^dydr + C fJ \{Q^y^,U^^,Q^y\U^^rdydr 

Kiei)G? 



-d^dydr + Ch,T 



(4.1.81) 



(4.1.82) 



(4.1.83) 



Moreover, Lemma 12.51 and Cauchy 's inequality and wave interaction estimates imply 

£ I /<•'' + ^-TifJI TI^**"^ + ^X' IK-^- 

+C6 [ \\{cb,^,0\\^dT + Ch,T e^- (4.1.84) 



Integrating (|4.1.80p with respect to ^, y and r and using p.23p . (|4.1.8ip - (|4.1.84p and the smallness of 
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X, 5, £ yield that 



<C5 J^^ \\{c^,ij,C}fdr + C J2 j[w'MX)fdT + C J !^\Gy\^d^dydT + Ch,T ei . 

e |q'| — 1 e e 

(4.1.85) 

On the other hand, from the fluid-type system p.6p . we can get an estimate for || (5'i-, VKr)|P as follows. 



^ \\i^r,Wr)fdT<cj^ \\{'S>y,^y,Wy,^yy,Cy)fdT + CS \\ ,/[U^\{^Sj , W)^ dT 

+Cs J^^ \\{^,W)fdT + Ci6 + CH^Tx) m J ^IGipd^dydr (4.1.86) 



[ [ l^ld'^'GMdydT + C.Tei 



\a'\ = l 

From (13.51) . we have 



||Cll'rfT<C ||VFy||2dr + C J\^y\UydT + C JlUy^ffdydr 
<C jlwWyW'dr + Cx' jl\\^y\?dr + C5 \\^\U^m'dT + Cu,t e m^dr. 



(4.1.87) 



In summary, collecting the estimates (|4.1.66p . (|4.1.79p . (|4.1.85p -( |4.1.87p . we complete the proof of Propo- 
sition 13.11 



4.2 Proof of Proposition 13.21 

Proof of Proposition 3.2. The proof is divided into the following five steps. 
Step 1. Estimation on || (</>, -0, C)('''> Oll^- 



V 2C 

Similar to (|4.1.10p . we muhiply ([535]) ^ by (j), ([X^ . ^V-i, dSSSa by V,, (1535)) ^ by respec- 

tively and adding them together to have 

{^^ a^^^l ^ C \ ^4/i(e) 2 ^Amg) 2 ^2«;(e) 2 
I Y + 2^"^! + 2^ T + 3z^ ] + ^r'^^y + 2^ + ^^^v 

\ i=2 / T i=2 

V - ' - 

= /6(0, V, C, '/'y, Cy) + ^V^l(A^5 - Qly ~ Qly) + MN^+4 " Q^y " Q^y) (4.2.1) 

4=2 

3 



[N8 - Qiy - Qiy + J2 ^'^^-V + ^^2/)] + ^3 + (' ' ' 



where 



2 4/i(e) t/ Aj(e) 2 7' 



(4.2.2) 
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and 



+^ [ - ^ / 6ieP(ni - nfi^ - nf'),de + E ^» / 66ni,de (4.2.3) 

i— 1 

3 „ 



Direct calculation shows that 



h <C\iVr,Zr,Vy,Uy, ey,Zy)\[Uy,,PyXv)\^ + , ^, C)|'] 

< C{5 + Ch,Te--my,^PyXy)\^ + C{5 + Ch^Te-- m,^, C)p. 



(4.2.4) 



Thus, integrating (|4.2.ip with respect to t and y and using Cauchy inequality yield that 

||(</.,V,C)(t,-)I|'+ r \\{i^yXy)fdT <C{5 + CH.Te-^) f Uvfdr 

+C{P + 5 + Cu.T £^) m^,OfdT + Ci, m Y.^Q% + Q%)dydT (4.2.5) 

By (P3T|) and ([XM]) . we have 

r / E(Q?y + Q',)rf2/rf^ < (4.2.6) 
From ([Qa]) . ([Oil) and ([OS]) . we have 

l(V',C)l[l^5|+El^'^+4l + l^8l]^y'^^^^''^^^ / ||(0,V,C,'/',>y,C,,V'w,C,,)ll'rfr. (4.2.7) 



Now we estimate the microscopic term J J K^dydr in (14.2.51) . We only estimate study K^i := 

~ L I ^"^^ / ^^^^^ ^ ~ Ilf')yd^dydT and X32 := ^ E / ^^Hij^d^ because the other 

terms J J K^dydr can be estimated similarly. 

E 

For K^i, integration by parts with respect to y and the Cauchy inequality yield 



K31 = j[ I {^My J e?(ni - nfi^ - Uf^)d^dydr 



UiyfdT + c{5 + Ch^Te"^) UifdT + Cp j\ y"6'(ni-nfi^-nf^)df|2d2/dT, 

(4.2.8) 



</3 

where the last term has been estimated in (|4.1.78p . 
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For we have 



K: 



32 



i=l 
3 



^^E [J {\CyM + \MC\ + \Zy 

■ [1 / ae.(ni - nfi^ - nf')de| + 1 1 ae.(nfi^ + nf^)dci 



dydr 



(4.2.9) 



Substituting (|4.2.6p - (j4.2.9p and (|4.1.78p into (|4.2.5p and choosing ;5,ir,(5, x suitably smah yield that 



ii(0,v,c)(r,-)ir+ \\{^jy,CyWdT<c{5 + CH,Tx) mvdT 

+Cp{6 + Ch^TX) / \\i^,^,CWdT + Ch,TX \\{(^y,^yy,CyyWdT + Ch,T,p 



|a'| = r 



(4.2.10) 



Step 2. Estimation on ||0j^(r, •)|| 



To estimate the term J \\<j)y\\'^dT, we firstly rewrite the equation p. 391) 9 as 



4/z(e) 



ii^)y^.y + i,Cy^H, 



e^Tll - nfi^ - nf')yd( -N5+ Qly + Qly, 



(4.2.11) 



by using the equation of conservation of the mass p. 39^ 1 . Multiplying ()4.2.1ip by we get 



(4.2.12) 



Integrating (|4.2.12p with respect to r, y and using (14.2.61) , (j4.2.7p and the Cauchy inequality yield 



Uyir,-W + UyfdT<C\\Mr,-)f + C \\ii^y,Cy)fdT 
+ C{6 + Ch^TX) jl\mi>,CWdT + Cx jlUlyyfdT 

+ Ch.t + f, J\J ^i(ni - nfi^ - uf^d^dydT. 



(4.2.13) 
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For the microscopic term J J \ J £,l{Ui - U'^f - Uf=')yd^\'^dydT, by (|4.1.42l) . we have 

e 

<C[J' J\J ei(L^Gr)yd^\^dydr + /J/l/ eiii^m^G,y)yd^\'dydT 



Then we have 

< c Y. nj^-:^\a'ofdidyLci^.c.,.)rjj^-^^did,d, 

\a\=2 e * * 

Similar estimates hold for K\ {i = 2,3). Moreover, 

Kt <C r f \ f e?Lj^'Q(Gi,Gi,)d^|2dydT 



+C j[ j\j ei^m{Q{'^mQ(Gi,G,)My)}d£.?dydT 



and 



<c 1^ [ [ \Ml±}:^Lt\:I^d^dydr 

<C{S + Ch,TX) I \\{(l),'>PX,<l^y,'iPyXy,<l^yy,'iPyyXyy)\\'^dT + Ch,T 



Substituting (|4.2.15p - (|4.2.18p into (|4.2.14p gives 
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(4.2.20) 



Thus we have 

Uv{r,-W + jlUyrdT<CU,{T,-)r + C jjUyXyWdr 

~ r ~ 1 

We now turn to the time derivatives. To estimate || (0t, ?At, Ci-)|P, we need to use the system p.27p . 
By muhiplying ((XTfl i by 0r, (|5^ 9 by -^ir, dSSZlla by tpir [i = 2, 3) and ^^TF^ a by Cr respectively, and 
adding them together, after integrating with respect to r and y, we have 



\\{<Pr,^rXr){T,-)rdT<C{5 + Ch,Te-^) / 1 1 (0, V, C) 1 1' + C / || (0„ ^„ C,) ir^T 

(4 22) 

+C m I !^\Gy\'dCdydT + Ch,Te^. 
Step 3. Estimation on \\{<pyj'4'yjCy){'''T')\\'^- 

Multiplying ^3M\i ^ by (j)y, ^3A6^ ^ by ^-ipiy, (EZSIa by and ^Mh ^ by gfrCi/, adding them 
together gives 

^2 +2Z^i^ ^ 2 3Z2^" 3Z '^12'!' + 2^ y V^yy + ^z^y'^yy 

T/" / ^ 

= l7((l),1p,C,(l}y,'^yXy,'^yyXyy) " (^5 - Qly - Qly){—^)y - ~ Q^V ~ Qw)'^^yy (4-2.22) 

1=2 

2r - ^ - 



'3Z2 
where 



l7{(l),ip,C,(f'y, i^y , Cy , V'yy , C; 



J/ J/-' 



and 



Y 1 2 V ^ 

~3§2^f^8 - {^^^^^^)y{^)yi^iv - (^l^)y^iy^V'iwy (4.2.23) 

_ V(^) ^. ^. _ i^^) ( J_) C - (^) C — C 
/ A lyViyViyy \ y )y\r.^2'y^y \ y /y'^J' 0^2 ^J'J' 

i=2 

< /3|(0yj,,Vaa,Cw)l^ + C/j((5 + C,,,T£3)|(<?!),V,C,0y,V'y,Cy)P, 



= {^)y I ei'(ni - nfi^ - nf'),de + I a6(ni - nf^^ - nf^),de 

(=2 

+(|%)i;[/ei^(ni-nfi^-nf^),de-E^» / ^i^^^i^de (4.2.24) 

i=l 

3 „ 

- E y ae.(ni - nfi^ - nf-)yd^] . 
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Integrating (|4.2.22l) with respect to T,y, and substituting (|4.2.6p . (|4.2.7|) and (|4. 2.191) into (|4.2.22p 
and (|4.2.24p . choosing f3,x,s smah enough, we have 



ll(</'y,V'2/,Ca)(T", Oil + \\{^yyXyy)\\ dr 

<P I UvyfdT + Cp{5 + Ch,Te^) I \\{(j),^pXAv.^vXy)fdT + Ch.Te-^ 



+CKTxjl j J '^IG.l'd^dydr + CY: fj J '^\d'^G\'d^dydr 



(4.2.25) 



|a|=2' 

+C,,TX fJJ ^\d"'G\'d^dydr. 

\a'\ = l e 



Again, to recover H^yylP in the dissipation rate, applying dy to (|3.27p 9. we get 

^lyr + iP~P)yy^-l[^Uly-^^^^U^y^]^^-J Cl^^Il?i''yydC-Qlyy-Qlyy^J CfGyyd^ (4.2.26) 



Note that 

iP - P)yy = ~z;^yy + ^Cyy ~ -{p " P)Vyy - ^Pyy --^{p~P)y~ -^<l>y. (4.2.27) 

(J O (J (J U (J u 

Muhiplying (j4.2.26p by -(j)yy and using (|4.2.27p imply 



i^iy<Pyy{T,y)dy + / / ^(/ylydydr 



<C{6 + Ch,TX) j^^ \\{<l>,i^X.(l)y,^yXy)fdT (4.2.28) 

+C jj\{i^iyyXyy)fdT + Cu,T e-^+cY. m J ^\d"G\'d^dydT, 



|a|=2 

where we have used the fact that 

|2 I I/O |2 



(IQlyyl + IQlyylndydT < Ch,T . (4.2.29) 

To estimate || (<^i/t, V'j/t, CyT)|P and || (</>,-,-, ^^..r) IP, we use the system p.27p again. Applying dy to 
p.27p . and multiplying the four equations of p.27p by 4>yr, 4'iyT, ipiyr (* = 2,3), Cj,t respectively, then 
adding them together and integrating with respect to r and y, we have 



|a|=2 

Similarly, we can obtain 



\\i4'yT,i'yT,CyT)\\ dr < C \\i4>yy,4'yy,Cyy)\\ dT + Ch,T£^ 

+Ch^TX [ \\{cf>y,^Jy,CyWdT + C{5 + Ch.Tei) [ ||(0,^,C)|pdT (4.2.30) 



+^ Hill '^\d'^^Mdyd. 



\\{(l)rT,lpTT,C,TT)\\ dr <C I \\{(j)yr,^yT,C,yr)\\ dr + Ch.T 

Jh. 

+C^,TX E /j|5"'(0,V,C)lPdT + C(<5 + C,,Te^) |j|(0,V,C)frf^ (4.2.31) 



|a'| = l' 



|a|=2 
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A suitable linear combination of (j4.2.25l) - (|4. 2.311) gives 



< 



\a\=2 
+ Ch,TX 



\Gifd^dydT + C{S + Ch^Te--) I ||(</), V,C)f 



(4.2.32) 



|a'|=l-^T 



C)fdT + C\T 



Step 4: Estimation on the non-fluid component: 



To close the a priori estimate, we also need to estimate the derivatives on the non-fluid component 
G, i.e., 9"G,(|a| = 1,2). For this, from (|3?TT|) . we obtain 

(4.2.33) 

Applying dy on (|4.2.33p . we have 

vGyr - («LmG), = |uiG„ - Pi(aGy) - i;[iPi(aM,) - -Lpf3(^iM^^) 

+2vQ{G,G'^'>) +vQiG,G) + vJi\ -VyGr- (4.2.34) 

Multiplying (|4.2.34p by ^ gives 



vG^ 



G,i 



2MW 



^(i;LmG), = ^ i^iG, - Pi(aG,) - V -Pi(eiM,) - —Pf^{^,M^^^) 



G, 



G, 



G? 



+2«Q(G, G^-^) + vQiG, G) + «Ji I ^ - ^'.G.^ + i;,^. 

Then Lemma l3.2l -Lemma 13.51 and Cauchy inequality give 

G ~ 
:^{vLmG) yd^dydr 

^ (vLMGy + 2vyQ{M, G) + 2vQ{My, G)) d^dydr 



> 



7a 



'^^^^^d^dydr - C{6 + C^.TX) 



d^dydr 



Ce / \Uy,^l^y,(:y)fdT-CH,Te--, 



and 



G„ 



Gl 



-^y 'Gr^ + + [uiGy - P,{^,Gy))y^d^dydT 

Kiei)p^2 



< 



32 



Kiei)G 



^d^dydT + C(,5 + C^,TX) 



Gid^dydT 



-C 



— — — dfdydr. 



(4.2.35) 



(4.2.36) 



(4.2.37) 
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Note that 

T ^ 



2v0 



20 



2V'^^eSs 1 LSI-- V 29^3 
Then Lemma 13.21 -Lemma 13.51 and wave interaction estimates imply that 

Q 

-^^^^-JLd^dydT + C{S + Ch,tx) 



(4.2.38) 



< 



and 



32 I h 



+C\6 + Ch,Tx) I \\{cb,ij,0\\m(dv)dT + Ch,Tei 



{vQ{G,G))^^J^dydT 



d^dydr 



(4.2.39) 



{vyQ{G, G) + 2z;Q(G, G,)) ^/Cdydr 



32 I h 



^ d^dydr + C{x + 5) 



d^dydr. (4.2.40) 



Thus, integrating (|4.2.35p with respect to ^, y and r and using (|4.2.36p . (|4.2.37p and (|4.2.39p . (|4.2.40p . 
we obtain 



^^''(T,y,Odedy + 



2M, 

+C{5 + C^,TX) 



Kiel) IP ,2 



G.yp(i,ed2/rf^ 



^|G,,pd^dydr + C ll('/'y„Cyy)ll'<iT 



Kiel) 



{\Gi\'' + \Gr?)d^dydT 



+C{5 + Ch,Tx) / [IK.^^,^!' + ll(0,V',C),ll']dT. 



(4.2.41) 



Similarly, 



IG. 



2M, 

< Ch.ts^ + C 



{T,y,^)d^dy + 

Kiel) 



Kiei),^ ,2 



+C(5 + Ch^TX) 



Kiel) 



Gyr\^d^dydT + C I \\{4>yr,CvrWdT 



(|Gi|^ + |Gj,|^)dedydr 



+ CC{6 + Ch,TX) i [||(0,^,C)lP + ll('^,V',C)rf]rfT. 
Step 5: Highest order estimates: 



(4.2.42) 
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Finally, we estimate the highest order derivatives, that is, — J ^piy't'yydy and 

/ / / M ^rfCc'yc^T" with |a| =2 in gX321)- To do so, it is sufficient to study / / m ^^^^ d'^ 

2) in view of p.38p and p.37p . Using the same idea in [19], we obtain the estimation for the highest 
order derivative terms, i.e. 



|a|— 2 2 ^ 

<C(r;o + <5 + Q,,Tx) I ^(|G,|2 + |G.|2 + |Gi|2Kdydr 

e 

2 

+Cirio + S + Ch,TX)Y. J, \\d"{<l>,i^,CWdT + Ch,T ei, (4.2.43) 



|a|=0' 

where 770 is defined in Lemma 13.31 
Noting that 

- j i^lytPyydy < ^HV-lylP + C'^ll^yylT 
and combining the estimates ()4.2.32p . (|4.2.4ip . (|4.2.42p and ()4.2.43p . we complete the proof of Proposition 
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